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Abstract 

A cubic algebraic equation for the effective parametrizations of the standard 
gravitational Lagrangian has been obtained without applying any variational prin- 
ciple. It was suggested that such an equation may find application in gravity theory, 
brane, string and Rundall-Sundrum theories. The obtained algebraic equation was 
brought by means of a linear-fractional transformation to a parametrizable form, 
expressed through the elliptic Weierstrass function, which was proved to satisfy the 
standard parametrizable form, but with g<i and g$ functions of a complex variable 
S^ instead of the definite complex numbers (known from the usual (arithmetic) the- 

ory of elliptic functions and curves) . The generally divergent (two) infinite sums of 
the inverse first and second powers of the poles in the complex plane were shown 
to be convergent in the investigated particular case, and the case of the infinite 
point of the linear-fractional transformation was investigated. Some relations were 
found, which ensure the parametrization of the cubic equation in its general form 
with the Weierstrass function. 
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I. INTRODUCTION 



The synthesis of algebraic geometry and physics has been known for a long time, beginning 
from chiral Potts model, algebraic Bethe ansatz (for a review of these aspects see l ) 
and ending up with orbifold models of string compactification. 2 In the context of string 
theories, the application of algebraic curves, related to Fermat's theorem, has also been 
pointed out. 3 

Concerning gravitational physics, which is an inherent constituent of any string, brane 
or ADS - theories, any applications of the theory of algebraic curves are almost absent. In 
this aspect perhaps one of the most serious attempts was undertaken in the recent paper 
of Kraniotis and Whitehouse. 4 Based on a suitably chosen metric of an inhomogeneous 
cosmological model and introducing a pair of complex vaiables, the authors have suc- 
ceeded to obtain a nonlinear partial differential equation for a function, entering the trial 
solution of the field equations. The most peculiar and important feature of the obtained 
equation is that it can be parametrized by the well-known Weierstrass function (for a 
classical introduction in the theory of elliptic and Weierstrass functions see Ref.5-7 ). This 
convenient representation enabled the authors to express important physical quantities 
such as the Hubble constant and the scale factor through the Weierstrass and the Jacobi 
theta functions. In fact, an analogy has been used with examples on the motion of a body 
in the field of an central force, depending on the inverse powers of the radial distance r. 
The cases of certain inverse powers of r, when the solution of the trajectory equation is 
expressed in terms of elliptic and Weierstrass functions, have been classified in details. 8 

Three important conclusions immediately follow from the paper of Kraniotis and 
Whitehouse, 4 and they provide an impetus towards further investigations. The first two 
conclusions are correctly noted by the authors themselves: 1. Other cases may exist, 
when solutions of nonlinear equations of General Relativity might be expressed in terms 
of Weierstrass or theta functions 9 , associated with Riemann surfaces. 2. The differen- 
tial equations of General Relativity in a much broader context might be related to the 
mathematical theory of elliptic curves and modular forms (for an introduction, see Refs. 
10-12 ) and even to the famous Taniyama-Shimura conjecture, stating that every elliptic 
curve over the field of rational numbers is a modular one. For a short review of some of 
the recent developments in the arithmetic theory of elliptic curves, the interested reader 
may consult also the monograph 13 . In fact, the eventual connection of General Relativ- 
ity Theory with Number and Elliptic Functions Theory was formulated even in the form 
of a conjecture that " all nonlinear exact solutions of General Relativity with a non-zero 
cosmological constant A can be given in terms of the Weierstrass Jacobi Modular Form" . 4 
Of course, such a conjecture is expressed for the first time, and yet there are no other 
solutions derived in terms of elliptic functions, not to speak about any classification of 
the solutions on that bases. 

The present paper will not have the purpose to present any new solutions of the 
Einstein's equations by applying elliptic functions, nor will give any new physical inter- 
pretation, which in principle should be grounded on previously developed mathematical 



techniques. Rather than that, in this paper an essential algebraic "feature" of the gravi- 
tational Lagrangian will be proved, which is inherent in its structure, mostly in its partial 
derivatives. This 'algebraically inherent structure" represents the third conclusion, which 
in a sense may be related to the problems, discussed in Ref. 4. 

However, this algebraic feature will become evident under some special assumptions. 
While in standard gravitational theory it is usually assumed that the metric tensor has 
an inverse one, in the so called theories of spaces with covariant and contravariant met- 
rics (and affine connections) 141 instead of an inverse metric tensor one may have another 
contravariant tensor g^ k , satisfying the relation gi$ k = l k 7^ 6 i . But then, since l k can- 
not be determined from any physical considerations and at the same time the important 
mathematical structure from a physical point of view is the Gravitational Lagrangian, a 
natural question arises: Is it possible that in such a theory with a more general 
assumption in respect to the contravariant metric tensor, the Gravitational 
Lagrangian is the same (scalar density) as in the usual case, provided also 
that the usual connection and the Ricci tensor are also given? Prom a phys- 
ical point of view, this is the central problem, treated in this paper, and the 
answer, which is given, is affirmative. Namely, it has been shown that if e, are 
the components of the covariant basic vectors, and dX^ are the components of a con- 
travariant vector field (which, however, are not contravariant basic vectors and therefore 
CidX = ^ 7^ 8\ ), then they satisfy a cubic algebraic equation. Of course, if dX % 
are to be found from this equation, then it can be shown that g^ k will also be known be- 
cause of the relation g^ k = dX^dX k . Also, it has been assumed that the affine connection 
T k j and the Ricci tensor Rij, determined in the standard way through the inverse metric 
tensor are known. 

The obtained cubic algebraic equation can be expressed in a very simple form, but 
unfortunately it is not easy at all to solve it. That is why a mathematical approach 
for dealing with such an equation has been developed, on which any further physical 
application will be based. The equation has been derived in two cases - when d 2 X l = 
and when d 2 X t 7^ 0. As will be shown, the first assumption means that dX has zero- 
vorticity components (and non-zero divergency, however), and in physical considerations 
this restriction can be imposed. The second assumption would mean that dX has both 
non-zero divergency and non-zero vorticity components. From the mathematical theory 
of cubic equations, the investigation of the two cases will not be different, because in the 
second case only the algebraic variety from the first case (with dX 1 ) will be supplemented 
with the components d 2 X l . It is worth mentioning also that the derivation of the equation 
does not presume zero-covariant derivatives of the covariant tensor field, thus leaving an 
opportunity to unvestigate the different kinds of transports on the space-time manifold. 
The algebraic equation may enable one to make a kind of a classification (from an algebraic 
point of view) of the contravariant tensors, satisfying the same gravitational Lagrangian. 

So far, the problem investigated here may seem to be of pure "mathematical" in- 
terest, but it may have also numerous physical applications. In supergravity theories, 
ADS/CFT, five-dimensional and brane theories 15,16,17 , one deals with an action, consist- 
ing of a gravitational part, added to a (for example) string action of the kind S str . = 



— -| J d 2 ^^/—hh al3 d a X^df3X^, where X^ are the string coordinates, h a/3 -the world sheet 
metric tensor, T- the string tension and the partial derivatives d a are taken in respect to 
the world sheet coordinates £ a = (t,ct). One can easily guess that the above described 
methodology can easily be applied to the string part of the action. More concretely, h a @ 
may be expressed as h a/3 = d£ a d^ , the gravitational metric tensor may be assumed to de- 
pend on the string coordinates and the derivatives of the string coordinates will be taken 
in respect to the world sheet coordinates C, a . Also, the partial derivatives in the gravita- 
tional part of the action may also be taken in respect to the coordinates C, a . As a result, 
taking the gravitational and the string part of the action together and without applying 
any variational principle, one would get the same kind of a cubic algebraic equation as the 
one, which will be proposed in this paper. In a sense, this dependence of g^ on the string 
coordinates is a sort of a coupling between the gravitational part of the action and the 
string one, and the resulting cubic equation may be called " an algebraic equation for the 
effective parametrization of the total Lagrangian in terms of the string coordinates" . The 
"coupling" between the two parts of the action provides another interesting possibility, if 
the first variation of the Lagrangian is performed, even without taking into account any 
equations of motion. Provided that the gravitational Lagrangian depends also on the first 
and second differentials of the metric tensor, the first variation of the Lagrangian can be 
regarded also as an cubic algebraic equation in respect to the differentials of the vector 
field 18 . 

There is also another interesting problem, which is related to the current trends in 
ADS/CFT correspondence 33 and the WZW theory of strings on a curved background. 34 ' 35 
For example, in WZW theory it is not clear how to relate the two-dimensional string world- 
sheet symmetries to the global symmetries (global coordinates) of the three-dimensional 
ADS spacetime, in terms of which the parametrization of the group element is pre- 
sented. If one has the ADS metric (and the ADS hyperboloid equation), then usually 
some parametrization of the global ADS coordinates is chosen, in terms of which the 
hyperboloid equation is satisfied. However, this is performed from the viewpoint of max- 
imum simplicity and convenience. The formalism, developed in this paper on the base 
of cubic algebraic curves gives a possibility to find some other parametrizations. For 
example, the (three-dimensional) ADS spacetime has a boundary which can be found 
when one of the (global) coordinates tends to infinity (for example r — > oo), and thus 
a two-dimensional space is obtained, which can be identified with the one, on which 
the world-sheet coordinates are defined. By "identification" one may mean that not the 
two-dimensional coordinates, but only the Ricci tensors and the Christoffel connections of 
the two space-times can be identified. Then, and as will further be shown in the general 
case, the obtained in this paper cubic algebraic equation will give an opportunity to relate 
one of the differentials with the Weierstrass function and the other differential - with its 
derivative. In such a way, one obtains a system of two equations in partial derivatives, 
from where the parametrizations can be found. In principle, 2+1 dimensional gravity and 
the WZW model of strings on an ADS background are very convenient for application 
of the algebro-geometric approach. Another interesting moment in these theories is that 
very often one has to deal with the ADS metric, written in different coordinates (includ- 



ing the two-dimensional coordinates of the worldsheet), and it may be supposed that the 
transition from one system of coordinates to another can probably be given by a linear- 
fractional transformation, which will be invesigated in this paper. From this algebraic 
point of view, it is interesting to investigate the coordinate transformations in Ref. 36 
and in Ref. 37. 

Until now, only the physical aspects of the implementation of the algebraic approach 
have been discussed. The mathematical theory of cubic algebraic equations is also worth 
mentioning, and creating the relevent mathematical methods will be the main purpose 
of this paper. In principle, the theory of cubic algebraic equations and surfaces has 
been an widely investigated subject for a long time. The mathematical theory of cu- 
bic hypersurfaces 19 puts the emphasis especially on the classification of points on the 
cubic hypersurface, minimal cubic surfaces, two-dimensional birational geometry and 
quasi-groups. But no concrete applications of cubic curves are given. In well-known 
monographs, 20 the general theory of afline and projective varieties and algebraic and pro- 
jective plane curves is exposed, and some examples are considered too, but the theory 
of cubic forms is restricted only with the Pascal's theorem. A more comprehensive in- 
troduction to the algebraic theory of second and third-rank curves, their normal forms, 
turning points (where the second derivatives of the curves's equation equal to zero), ra- 
tional transformations and etc. is given in the book of Walker 21 . Of particular relevence 
to the present research will be the theorem 21 that if f(x, y) = is a non-degenerate cubic 
curve, then by introducing an afnne set of coordinates xi = -, y± = - and choosing the 
turning point at (0, 0, 1), the curve can be brought to the form y 2 = g(x), where g(x) is a 
third-rank polynomial with different roots. However, the situation is much more interest- 
ing in the complex plane, where one may define the lattice A = {muj\ + nuj2 \ m,n e Z; 
u>i, 0J2 E C, ira— > 0}. Let a mapping / : C/A — ► CP 2 is defined of the factorized along 
the points of the lattice part of the complex plane into the two-dimensional complex 
projective space CP 2 . If under this mapping the complex coordinates z are mapped as 
z — > (p(z),p (z), 1) when z ^0 and z — * (0, 1,0) when z = (p(z) is the Weierstrass 
elliptic function), then the mapping / maps the torus C/A into the following affine 
curve y 2 = 4x 3 — g^x — #3, where gi and g^ are complex numbers. 13 The important mean- 
ing of this statement is that excluding the points on the lattice which may be mapped 
into one point of the torus (where the Weierstrass function has real values), the mapping 
z — > (x, y) = (p(z), p (z)) parametrizes the cubic curve. The consequence from that is also 
essential since the solution of the resulting differential equation can be obtained in terms 
of elliptic functions. 5,13 ' 22 In spite of the fact that the parametrization can be presented 
in a purely algebraic manner, it is inherently connected to basic notions from algebraic 
geometry such as divisors and the Riemann-Roch theorem, 23 which reveals the dimension 
of the vector space of meromorphic functions, having a pole of order at most n at the point 
z = 0. This should be kept in mind because some results may be obtained by algebraic 
methods only, but the explanation may probably be found in algebraic geometry. In the 
present paper, a more general parametrization of a cubic curve is considered , when g% and 
#3 are not complex numbers (the so called Eisenstein series g<i = 60 ^ ^4 = ^ , ^ r )4 i 

wcr n,m 



#3 = 140 ^2 ^g — Yl (ra+mri 6 )' but com P lex functions. It has been proved that if 

wcr n,m 

the Weierstrass function parametrizes again the cubic curve, then the infinite 

sums (in pole number terms) ^ -^ for n = 1 and n = 2 turn out to be finite 

wcr 
(convergent) ones, in spite of the fact that in the general case they might be infinite 

ones (divergent). The explanation of this fact from the point of view of algebraic geom- 
etry remains an open problem, but it can be supposed that standard arithmetical theory 
of elliptic functions and algebraic equations is contained in some other, more general 
theory, which may be called non- arithmetical theory, and from this theory the standard 
parametrization should also follow. The considered case of parametrization of a cubic 
curve with coefficient functions of a complex variable, although performed in this paper 
in a trivial algebraic manner, is the first step towards constructing such a theory. At 
least, a certain motivation from a physical point of view is evident for constructing such 
a theory. 

The above-presented outlook on standard parametrization implied the use of affinc 
coordinates, which unfortunately exclude from consideration the infinity point. But the 
infinity point cannot be ruled out not only from mathematical grounds, but also from phys- 
ical considerations. For example, in the five- dimensional Randall-Sundrum model, 29 ' 30 
one has to assume a compactification into a four- dimensional universe from an infinite 
extra dimension, containing also the infinity point. From this point of view, the more 
convenient transformation, chosen in the present paper, which brings the cubic curve into 
a parametrizable form, is the linear-fractional transformation. This transformation allows 
one to parametrize with the Weierstrass function the ratio of the two of the parameters, 
entering the linear-fractional transformation and in the case the parameters in this trans- 
formation represent complex functions. Of course, the other parameters remain unfixed, 
leaving the opportunity to determine them in an appropriate way. In a sense, from most 
general grounds the appearence of the Weierstrass function in the linear-fractional trans- 
formation might be expected , since according to a theorem in the well-known monograph 
of Courant and Hurwitz, 24 an algebraic curve of the kind w 2 = a v 4 + aiv 3 + a2V 2 + a\v + a 4 
can be parametrized as v = a J'\ z J )+d = <p(z) and w = (p (z) by means of the transformations 
v = f^4^ and w = w\ (e " d ~^a ■ In the case clq = (which is the present case of a third-rank 
polynomial), tp(z) will be a linear function of p(z). 

In the present paper, however, the situation is quite different, since the linear-fractional 
transformation is applied only in respect to one of the variables (v), and in order to get the 
standard parametrizable form w 2 = 4v 3 — g2V — g% (with g 2 and g%— complex functions), 
an additional quadratic algebraic equation has to be satisfied. What is more interesting is 
that after the parametrization is performed, the linear-fractional transformation turns out 

l r ii-i M z )p{ z )+-+ B { z )-4;(p 2 { z )) i A n n n r r 

to be ol a more general kind v = —. . , c , . —. f d , where A, B, G, D are tunctions or 

to C(z)p (z)+D(z)+f ' ' ' ' 

2, and the expression for v represents a rational transformation of the kind v(z) = -^f\- 
Now from another point of view it can also be understood why it is justifiable to apply 
the rational transformation only in respect to v and not in respect to w. The reason is in 
a well-known theorem 13 from algebraic geometry that "each non-degenerate cubic curve 



does not admit a rational parametrization" . Since each non-degenerate cubic curve can 
be brought to the form w 2 = v(v — l)(v — A), (A 7^ 0, l),the essence of the above theorem 
is that this (algebraic) form cannot be satisfied by a rational parametrization of both 
v = t44 and w = -J4. 

Q(z) R(z) 

In respect to the problem about finiteness of the infinite sums ^ \ for n = 1 and 
n = 2, the application of the linear-fractional transformation has also turned out to be 
useful, when the case of poles at infinity is considered. In principle, in this paper two 
separate cases are distinguished - the first case of an infinite point of the linear-fractional 
transformation and the second case of poles at infinity, when in the infinite limit u — » 00 
the sum ^ -^ tends to the Riemann zeta function £(n). For this partial case and applying 
a different mathematical method, based on the Tauber's theorem, 34 one comes again to 
the fact about the finiteness of G\, proved in the general case by performing a Loran 
function decomposition. 

The present paper is organized as follows: 

Sect. II gives some basic formulaes about the so called gravitational theory with con- 
travariant and covariant metric tensors. In Sect. Ill the third-rank algebraic equation 
has been derived, starting from the standard gravitational Lagrangian. Also, the effective 
parametrization problem has been formulated in an algebraic language. In Sect. IV the 
general mathematical set-up for parametrization of the cubic equation has been discussed, 
and some physical motivation for the application of the linear-fractional transformation 
from the point of view of Randall-Sundrum theory has been presented. Sect. V shows how 
the cubic algebraic equation transforms under the action of the linear-fractional transfor- 
mation. Sect. VI shows how from the transformed cubic equation one can get the standard 
parametrizable form of the cubic equation (with gi and g% -complex numbers) and also 
the quadratic algebraic equation is derived, which has to be fulfilled if the parametrizable 
form holds. The approach is valid also when g 2 (z) and g%{z) are complex functions. In 
Sect. VII it was proved that the nonlinear and nonpolynomial transformation from the 
"unbar" to the "bar" variables is also invertible, thus giving the opportunity to write 
down two of the additionally imposed equations in terms of the new "bar" variables. In 
Sect. VIII the Loran's decomposition has been performed of the functions on the both 
sides of the algebraic equation (-£) 2 = M (z) p 3 + N (z) p 2 + P(z) p + E(z) , where p(z) is the 
Weierstrass function and M, N, P, E are functions of the complex variable z. A system 
of three iterative (depending on n) algebraic equations has been obtained, representing a 
necessary (but not sufficient!) condition for parametrization of a cubic equation of a gen- 
eral form with the Weierstrass function. It is not occasional that the condition is called ' 
a necessary, but not sufficient one" , because in principle more algebraic equations have to 
be solved in order to prove the existence of such a parametrization. In Sect. IX the possi- 
ble parametrization of the more simplified cubic equation [p (z)~\ = 4p 3 — gi{z)p — Qz{z) 
has been considered, and of course the main motivation for considering such a case is 
the close analogy with the well-known case, when g-i and g% are complex numbers. By 
calculating the coefficients in the negative power Loran expansion and combining them, 
it has been proved that the sums Yl ~ an d Y2 ^2 represent finite (convergent) quantities. 



The other equations for the other values of m — —3, —1 have been presented in Appendix 
A; for values of m = Ik - in Appendix B and for m = 2k + 1 and m = —k - in Appendix 
C. The equations in these Appendixes in fact complete the proof that all the Loran co- 
efficient functions can be uniquely expressed through a combination of the finite sums 
G n . The calculations are purely technical but they serve as a strict mathematical moti- 
vation and a proof of the new and basic qualitative fact that the Weierstrass function can 
parametrize the simplified form of the cubic equation with coefficient functions gi{z) and 
gz{z). This fact probably might represent one of the starting points in the creation of the 
so called non- arithmetical theory of algebraic equations. Sect.X investigates the positive- 
power decomposition of the above-mentioned equation, and from the convergency radius 
of the infinite sum the asymptotic behaviour of some of the Loran coefficients was found 
to be — j~i- Sect. XI considers the case of poles at infinity in the positive-power Loran 
decomposition, and from the requirement to have a certain convergency radius, expres- 
sions for some of the Loran coefficient functions are obtained. In Sect. XII a split of 
the original cubic equation into two equations is performed, and based on the fact from 
Sect. IX, it has been proved that the parametrization of the first equation leads to a 
parametrization of the second equation. For the two " splitted-up" equations, Sect. XIII 
presents an algebraic equation, defined on a Riemann surface, which has to be satisfied 
if the so called j— invariants of the two equations are to be equal. In Sect. XIV on the 
base of the Loran function decomposition of g2(z) an infinite sum is obtained in which the 
coefficient functions contain the sums G n . In Sect. XV this formulae has been combined 
with a proof that the Tauber's theorem can be applied, and this combination resulted in 
an expression for G\ in the limit of poles at infinity. In this partial case, the expression 
again confirms that the sum dis convergent. In Sect. XVI the case of the infinite point of 
the linear-fractional transformation is considered, and the approach essentially represents 
a combination of the "split-up" approach from Sect. XII and the approach from Sect. 
VI , based on the derivation of the additional quadratic equation. In Sect. XVII the 
relation between the two integration constants is found which appear in the process of 
integration of the two splitted-up equations. A peculiarity of the developed approach is 
the appropriate "fixing" of some of the functions in the linear-fractional transformation 
so that the simplest and most trivial form of the quadratic equation from Sect. VI is 
obtained. Sect. XVIII starts with an algebraic equation of a fourth rank, derived from 
the original equation in the case of an infinite point of the linear-fractional transforma- 
tion. The main result here is that the constant Weierstrass function can parametrize this 
equation if it should be fulfilled in the entire complex plane. For the same equation, Sect. 
XIX investigates the second case, when the fourth-rank algebraic equation determines a 
Riemann surface for the pair of variables (p(z) = W\(z) + iw^iz); z), and six values of w\ 
are found, satisfying this equation. Sect. XX finds the necessary and sufficient condition 
for parametrization with a constant Weierstrass function, based again on the approach 
of Riemann surfaces. As a result, an integrable nonlinear equation is obtained for the 
coefficient functions of the algebraic equation. The coefficient functions in the solution of 
this equation appear in powers of non-integer (fractional) numbers. 



II. Covariant and Contravariant Metric Tensor 



Usually in gravitational theory it is assumed that a local coordinate system can be defined 
so that to each metric tensor g^ an inverse one g^ k can be defined 

9ij g jk = 5 k = {0 if i^ k and 1 if i = k} . (1) 

However, the notion of a reference frame can be defined in different ways in Ref. 25 - 
coordinate, tetrad and monad. In the last case the contravariant vector field dx % of an 
observer, moving along a space-time trajectory, represents a reference system. In such a 
case one may have instead of (2) 

e i d^ = fl^S{ = S(e i ,<y). (2) 

In the context of the so called dual algebraic spaces in Ref. 26, S(ei, dx j ) is called a 
contraction operator. Assuming that an inverse operator of contraction /■ exists, it can 
easily be obtained, as in Ref. 14 

e j = fidx 1 . (3) 

Therefore, the metric tensor field g can be decomposed in respect to the conravariant 
basic eigenvectors in the following way 

9 = ftj(e* <g> e j ) = g ij f k f l dx k dx\e k <g> ej) = (dx k dx l )(e k <g> ej), (4) 

and the contravariant components g^ of the tensor field g are represented as a contraction 
of the two vector fields dx l and dx^\ 

g lj = dx l dx 3 . (5) 

It is important to realize that this definition of a contravariant tensor field is not related 
to any notion of infinitesimality. In order to understand this, consider a set of global 
coordinates X M , defined on the given manifold and depending also on some other (local) 
coordinates. Then the set of global coordinates, regarded as functions of the local ones, 
can be considered as a system of equations, defining some algebraic surface. Provided 
that the partial derivatives of the global coordinates in respect to the local ones are non- 
zero, at each point of this surface the corresponding tangent space can be determined, 
and the differentials of the global coordinates are defined on this tangent space. If one 
assumes that the global differentials are infinitesimally small, then either the (partial) 
derivatives of the global coordinates or the "local' differentials should be small. However, 
the parial derivatives cannot be small, because one considers arbitrary global and local 
coordinates on the manifold. Also, if the local differentials are assumed to be small, then 
they will not be allowed to take arbitrary values. But this will mean that a large variety 
of integral curves on the manifold should be excluded from consideration. This will be 
unacceptable since one would like to define integral curves through each point of the 



manifold and moreover, it would contradict to our initial assumption about existence of 
a tangent space at each point of the surface (or manifold). Therefore, as a partial case, 
each of the local differentials should be allowed to take arbitrary numerical values and of 
course, they may be equal also to an arbitrary function of the local coordinates. 

It might be concluded, therefore, that since the partial derivatives and the local dif- 
ferentials cannot be infinitesimally small, then the global differentials cannot be also 
infinitesimally small. 

Apart from the definition (5) of a contravariant tensor field, we have also the definition 
of a length interval in Riemannian geometry 

ds 2 = l 2 (r) = g ij dx i dx j (6) 

If we would like to "incorporate" in this definition the standard definition of an inverse 
metric tensor as gijg^ k = 8^ , we can set up for the ordinary inverse metric tensor 

g* = ^dx l dx j . (7) 

V 

Therefore, in terms of the differentials, the ordinary inverse tensor g^ can be represented 
in the same way as in (5), but divided by the length interval. However, usually the lenght 
interval is not known, so from a physical point of view the definition (7) is undesirable 
and this is the motivation to deal further with the definition (5) of a contavariant tensor 
field. In order to distinguish the "newly" defined tensor in (5), a "tilda" sign has been 
placed. 

From (3) and (7) it follows 



2 ~~ 9klfi fj 



dx l dx j = 0. 



Clearly, the requirement for existence of an inverse contraction operator is equivalent to 
putting / = 1, i.e. assuming that there is a unit lenght interval, which is again physically 
unacceptable, and it is more natural to assume that the lenght interval is varying. Let 
us assume that l 2 and f\. are known in advance, then it can be be investigated which is 
the algebraic variety of values of dx l , satisfying this quadratic form. The main difficulty 
in this approach is that f l k cannot be determined from physical considerations. That is 
why the aim in the next section wil be to derive an algebraic equation, in which known 
physical quantities will enter - the metric tensor g^, the Christoffel connection Tf- and 
the Ricci tensor R^j. 
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III. CUBIC ALGEBRAIC EQUATION 
NOT FOLLOWING FROM 
A VARIATIONAL PRINCIPLE 

Further in this paper it shall be assumed that if X 1 are some generalized coordinates, de- 
fined on a n— dimensional manifold with coordinates on it (x 1 , x 2 , , x n ), then the differ- 
ential dX l is defined in the corresponding tangent space Tx of the generalized coordinates 
X 1 = X l (x 1 ,x 2 ,x 3 ...,x n ). Even if written with a small letter, it shall be understood that 
x 1 represent generalized coordinates. 

Our starting point for the derivation of the cubic equation will be the assumption that 
in spite of the choice for the contravariant metric tensor, the gravitational 
Lagrangian L = —^/—gR should be the same, provided also that the Ricci 
tensor does not change under the definition of the contravariant metric tensor. 
The meaning of this statement is the following. 

Essentially, the gravitational Lagrangian will have two representations. The first 
representation is based on the standardly defined Christoffell connection Y\ k 

T ik = 7;9 ls (9ks,i + 9is,k ~ 9ik,s) (9) 



T 



and the Ricci tensor 



fiV 1 fiV 1 

p ik il I -pi -pm pm-p! /i n\ 

8x l dx k L H L km ■ l iU J 

The second representation of the gravitational Lagrangian will be based on the def- 
inition (5) of the contravariant metric tensor cp fc = dx^dx k . Therefore, the Christoffell 
connection and the Ricci tensor will be different from the previous ones and will be de- 
noted respectively by Y\ k and R ik 

r^ = -j9 ls {9ks,i + 9is,k ~ 9ik,s) = -dx l dx s g ks>i + -dx l dx s g iS:k - -dx l dg ik , (11) 

~ fir 1 fir 1 ~ ~ ~ ~ 

p ik _ il _|_ -pi -pm -pm-pl ('{^ 

fix 1 fix k ik lm il km ' ^ ' 

The gravitational Lagrangian in this second representation is 
L 2 = -\/ ZI gR = -\f-gg xk Rik = -\f ZI gdx l dx 



fir 1 fir 1 

ikr> __ I -J^ij^k? UL ik UL il' 



fix 1 fix k 

— y—gdxdx (T ik T lm — T u T km ) . (13) 

Note that physical meaning of this Lagrangian will depend not only on the properties 
of the (covariant) metric tensor, but also on the first and the second differentials dx l and 
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d 2 x l . It should be mentioned also that the notion of a metric tensor, depending on gen- 
eralized coordinates, understood in the sense of a hypersurface (an infinite-dimensional 
manifold of all space-like hyper surf aces, embedded in a given Riemannian spacetime), has 
been introduced long time ako by Kuchar in Ref.31. In such an approach, the description 
of the gravitational field essentially depends on the tangential and normal deformations 
of the embedded hypersurface. In our case, we do not restrict to space-like hypersur- 
faces, but the notion of the differentials begins to play a self-consistent role, similarly to 
the dynamics and the deformations of the hypersurface in Kuchar 's approach. Yet, the 
standard gravitational physics with the usual inverse metric tensor is contained in the 
proposed in this paper approach, because one can identify the components of the usually 
known inverse metric tensor with the components of the contravariant metric tensor, de- 
fined in terms of the differentials. Thus one can obtain a system of first -order nonlinear 
differential equations in partial derivatives. The solution of this system may enable one 
to chose such global (generalized) coordinates, in terms of which the usual inverse tensor 
will be equivalent to the contravariant one in terms of the differentials. 

Let us now use expressions (5) for the contravariant metric tensor g^ and (11) for the 
Christoffel connection T^- in order to rewrite the gravitational Lagrangian in the second 
representation. The first two terms in (13) can be calculated to be 

-yf=gdx i dx k {— j f - — £) = ^f^gdx i dx k dx l {g is j \ X . } - -pg ik ,i + -gu„ '' "' ' - 



ox 1 ox" ox" 2 2 ox K 

= -^/^gdx i dx l {pT r il g kr dx k - T r ik gi r d 2 x k - T\ (i g k)r d 2 x k } , (14) 

where p is the scalar quantity 

,. d(dx l ) , . 

p = dtvidxj = — - . , (15) 

ox 1 

which "measures" the divergency of the vector field dx. It will be more interesting to 

calculate the contribution of the second term in (12) 

- v ^gdx i dx k (f l ik f]^ n - TilTkm) = -^V^gdx i dx k dx l dx m (-dgi m dx s g ks4 - dg ik dx r g mr ^+ 



+dgudx r g mr)k + dg km dx s gi S)i ) - y/^gdx l dx k dx l dx m dx s dx r [g ks ^g mr ^ - gi s ,ig mr , k ] =0 (16) 

and the first differential dg^ is represented as dg^ = -^dx s = T^gjydx 8 and T r si is the 
standard Christoffell connection. Therefore, the second two terms in (13) give no contri- 
bution to the gravitational Lagrangian. This is not surprising, since the "factorization" 
of the contravariant metric tensor as dx % dx^ introduces an additional "symmetry", due to 
which all the terms in (16) cancel. That is why the second representation of the grav- 
itational Lagrangian will be given only by the first two terms —^/—gdx l dx k {-Q^- — ^f ) 
in expression (13). 

12 



Concerning the first representation of the gravitational Lagrangian, it was based 
on the standard Christoffell connection T k - , the Ricci tensor R ik and the usual inverse 
metric tensor g l K The basic assumption at the beginning concerned the gravitational 
Lagrangian and the Ricci tensor, which means that together with the inverse metric tensor 
g % i , another contravariant tensor g^ = dx l dx^ exists, which enters the expression for 
the first representation of the gravitational Lagrangian 

Li = -^/=gg ik R ik = -V=gdx i dx k R ik . (17) 

Comparing this representation with the second one, given by expression (13) 



L2 = -V^gg 1 Ru = ~V^gdx l dx {pT r u g kr dx - T r ik g lr d x - Y r l{i g k)r d x } (18) 

and remembering the initial assumption, acccording to which the Lagrangian should 
be one and the same in both the representations (i.e. L\ = L2), one arrives at 
the following algebraic equation in respect to the first differential dx k and the second 
differential d 2 x k 

dx l dx l (pT r u g kr dx k - T r ik g lr d 2 x k - T r l{i g k ) r d 2 x k ) - dx i dx l R a = . (19) 

In the limit d 2 x k = this equation assumes the form of a manifestly cubic in respect 
to dx l algebraic equation 

dx l dx j dx k pT r j ( i g k y - R ij dx l dx j = 0. (20) 

Equation (20) is the basic equation, which shall be investigated further in this paper. 
Most importantly, it is manifestly cubic in the differentials dx % . Due to this reason, one 
qualitative argument can be given in favour of such a Lagrangian. In 1988, Witten derived 
the Lagrangian for 2+1 dimensional gravity in Ref. 32, which is also manifestly cubic in 
the chosen gauge variables A M . The Lagrangian was obtained under the assumption that 
there is an isomorphism between an abstractly introduced (d— dimensional) vector bundle 
with a structure group SO(d — l,d) and the tangent bundle of the given manifold, on 
which the metric is the induced one from the metric on the vector bundle. Besides, the 
verbein was assumed to be invertible, but as Witten remarks "permitting the verbein to 
not be invertible seems like a minor change." In the present case, we don't have at all any 
symmetry on the tangent bundle, neither is anything supposed about the dimensionality of 
spacetime or even about the existence of the usual inverse tensor, but yet the Lagrangian 
exibits the same cubic structure. Therefore, it may be concluded that the cubic structure 
of Chern-Simons theory 32 is inherent in the structure of the gravitational Lagrangian 
itself, and not in the additional assumptions in Ref. 32, which affect the choice of the 
gauge variables. In view of this, it might be interesting to investigate whether there is a 
transition from the Lagrangian in our case to the Lagrangian for 2+1 dimensional gravity, 
presented in Ref. 32. 

Of course, one might slightly modify the basic assumption, concerning the first rep- 
resentation of the gravitational Lagrangian. For example, instead of assuming that the 
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Ricci tensor will be the same in both representations, one might instead assume that the 
Ricci tensor should not change. In such a case again in the limit d 2 x k = the cubic 
algebraic equation will be in a form without the quadratic in dx l term 

dx i dx j dx k pT r j{i g k)r - R = . (21) 

One can write down also the vacuum Einstein's equations when the contravariant tensor 
is defined as g 1 ^ = dx l dx^ 

= --pgi j T r mn g kr dx k dx m dx n + -gij{Y r km g nr + T r n{m g k)r )d 2 x k dx m dx n + 

+ P T r ljgkr dx k - (T r ik9jr + T r j{l g k)r )d 2 x k . (22) 

Note the following subtle moment : since we have an expression equal to zero, this time 
it is not necesssary to assume that the above algebraic equation is valid under the 
assumption that the Ricci tensor does not change. Therefore, equation (22) provides the 
interesting possibility for classification of all solutions of the vacuum Einsteins equations 
with a given metric tensor gij and unknown contravariant tensor g iJ = dx % dxK In spite 
of the presence of the second differentials d 2 x k , equation (22) can be treated on an equal 
footing as an algebraic equation simply by "extending" the algebraic variety for the {dx k } 
variables with the new variable dy k = d 2 x k . However, if additionally it is assumed that 
the Ricci tensor does not change under the definition of the contravariant tensor (i.e. 
Rij = Rij), then one has 

(T r ik g jr + T r Kl9k)r ) d 2 x k = V Y\ 3 g kr dx k - Rij (23) 

and consequently all the terms with d 2 x k in the Einstein's vacuum equations (22) drop out 
and the algebraic equation becomes a cubic one in respect to the variables dx k only. The 
above analyses has the purpose to demonstrate that depending on the initial assumptions 
about the Ricci tensor or scalar curvature, the structure of the algebraic equation also 
changes. 

In an algebraic language, 20,27,28 the investigated problem can be formulated in the 
following way: 

Proposition 1 Let the differentials dx l (i = 1, ...., n; n is the space-time dimension) repre- 
sent elements of an algebraic variety X = (dx 1 ,dx 2 ,dx n ). For different metric tensors 

(and therefore - different connections Y k - and Riemannian tensors R- lk ), a set of polyno- 
mials (cubic algebraic equations) F(X) = may be obtained, which are defined on the 
algebraic variety X and belong to the ring R[dx x , dx 2 , ..., dx n ] of all third-rank polynomials. 
Then finding all the possible parametrizations of some introduced generalized coordinates 
X 1 (x 1 ,x 2 ,x 3 , ..,x n ) is equivalent to: 1. Finding all the elements dX l of the algebraic 
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variety X, satisfying the equation F(X) = 0. These elements will be represented in the 
following way 

(Jjilj ar I Jb - • • . Jb - (y^ 4 \ Jb - Jb . • • Jj I * J- :„"l Jb - Jb • • • • Jb I . -ilj£ -i I iX 1 - vC 1 . • . Jb ) ) ' \ "^ / 

2. Finding all the solutions of the above system of partial differential equations. 



In the present case, the algebraic equation is obtained before performing the variation 
of the Lagrangian, unlike the considered in Ref.18 another case, when again a cubic 
algebraic equation had been obtained after performing a variation. 

Let us comment briefly on the important from a physical point of view assumption 
d 2 x % = 0, under which the cubic equation (20) was derived. Suppose that for the set of 
generalized coordinates X 1 = X t (x 1 ,x 2 ...x n ) one has 

dX = aidx 1 (25) 

and let us assume that the Poincare's theorem is fulfilled in respect to dx l , i.e. d 2 x % = 0. 
Then 

d 2 X = d ai dx' + a t d 2 x l = ^dx j A dx l = (^ - ^Pr)dx i dx j . (26) 

OX 3 OX 1 ox 1 

Clearly, d 2 X = only in the following two cases: 1. <2j = const., i.e. dX 1 is a full 
differential. 

2. (rota)j-j = I s : — tA- = 0. The last means that if dx % are considered to be basic 

v i »j ox] ax% 

eigenvectors, then dX % have zero-vorticity components. Throughout the whole paper dX 1 
shall be considered as vector field's components in the tangent space Tx- 

Note also that the algebraic equation (19) with the first and the second differentials 
dx l and d 2 x % takes into account two important physical characteristics of the vector field 
dx l - the divergency p and the vorticity (through the term d 2 x l ). It might be required that 
these characteristics vanish, i.e. p = d 2 x l = 0. In such a case one is left only with the 
equation 

R ik dx i dx k = 0. (27) 

If additionally the requirement for the existence of the (usual) inverse metric tensor is 
imposed then the intersection variety of the quadratic form (27) with the quadratic forms 
(one-when Sj = 0,and the other - when 8\ — 1). 

g ik dx k dx j = SI . (28) 

has to be found. From the two last equations one easily obtains 

(R ik - \g %k R)dx k dx 3 = ~RS{ , (29) 
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in which the left-hand side is identically zero for every dx l in view of the Einstein' equations 
Rik — \g%kR = 0, but the right-hand side is zero only for i ^ j, but not also when % = j. 
Therefore, the Einstein's equations are obtained only in one case and not in the other 
case. In fact, it shouldn't be surprising that the Einstein's equations cannot be obtained 
for both the cases i ^ j and i = j. One should remember that the usual variational 
procedure in general relativity takes into account also the variation of the volume factor 
y/—g, while in our purely algebraic treatment and without any variation this volume 
factor was not subjected to any changes at all. Moreover, it is one standard procedure 
to perform the variational procedure with the usual gravitational Lagrangian and the 
inverse metric tensor (when the Einstein's equations are obtained) and its quite a different 
procedure to start from the other representation of the gravitational Lagrangian (where 
the variables to be variated are g^-, T^. (or g^k) and dx % and d 2 x l ) and afterwards to impose 
the requirement for identification of the contravariant metric tensor with the inverse one in 
the form of another, additional equation. So one should not even hope to obtain anything 
similar to the Einstein's equations. However, as already shown, if one has the Einstein's 
equations, one may still ask the question are they satisfied under the new definition of 
the contravariant tensor. 



IV. A GENERAL MATHEMATICAL SETUP 
FOR TREATING THE CUBIC 
ALGEBRAIC EQUATION (20) 

The subsequent investigation of equation (20) will be restricted to the case of a 5- 
dimensional space-time, although the approach of course may be applicable to any di- 
mensions. The main reason for chosing a 5-d spacetime is related to the widely discussed 
Randall-Sundrum (R - S) model, 29 ' 30 in which the process of compactification of the five- 
dimensional universe to our present four-dimensional universe is related to the existence 
of a large extra dimension. In the original R - S scenario the metric was chosen to be 

ds 2 = e- 2krcrr '^ v dx^dx v + r 2 c dxl (30) 

where r c is a compactification radius, rf v is the ordinary Minkowski metric, x$ C [0, 7r] 
is a periodic coordinate, \xv are four dimensional indices and A; is a scale of order of the 
Planck scale. Instead of the coordinate x$, one may chose for example a fifth coordinate 
X 5 = kr c x$, which in view of the largeness of the scale factor k may be assumed to range 
to infinity. But the infinity point, from a purely mathematical point of view, may be treated 
on an equal footing with all other points in the framework of projective geometry. 7 ' 21 ' 38 ' 39 
In the present case the infinity point shall be realized in respect to dx 5 after performing 
the linear-fractional transformation 

~5 

, 5 adx +b 
dx 5 = -^ , (31) 

cdx + d 
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where a, b , c and d will be assumed to be functions, depending on the complex variable 
z (or on two complex variables). Also, the remaining four-dimensional space-time with 
coordinates (x 1 , x 2 , x 3 , x 4 ) may be complexified in the following way 

Z1 — X1 + ix 2 ; z 2 = x 3 + ix 4 . (32) 

It is easily seen that the infinity point in respect to dx 5 is situated at dx 5 = — - and 
it is a zero point for the complex plane dx 5 . The convenience of the linear- fractional 
transformation from a physical point of view matches also the mathematical requirements 
of the problem. In order to parametrize the third-rank algebraic equation (15), written 
in a two-dimensional form, one has to bring it to the form 

{dx 5 ) 2 = 4{dx 4 ) 3 - g 2 {dx 4 ) - g 3 . (33) 

In the case when g 2 and g 3 are complex numbers g 2 = 6OG4 = 60 Yl -^ and g 3 = 140Gq = 
Y ^e ) standard algebraic geometry contains a well-known prescription how to parametrize 

U! 

this algebraic curve 13 by introducing the variables 

dx 4 = p(z) dx 5 =p'(z) , (34) 

where 2; is a complex variable and p(z) is a complex meromorphic function - the 
Weierstrass function 



p( z ) = h+H 



1 1 



z 
and the summation is over all non-null elements 



z — zu) 2 zu 2 



(35) 



w C A = {{mvui + nzu 2 ) \ m,n C Z (integer numbers), vo\ , va 2 C C, Im > 0} . 

Since further in the text the parametrization (34) will be repeatedly used, it is instructive 
to give just an idea how in classical textbooks it is proven that the parametrization (34) 
satisfies eq.(33). Let us take for example the proof, given in Ref.13, where the basic 
idea is to compare the Loran expansions for the non-positive degrees of z for the function 
[p (z)\ and for the polynomial ap 3 (z) +bp 2 (z) + cp(z) +d, where a, b, c and d are complex 
numbers. If the corresponding coefficients in the Loran expansion of these two expressions 
are equal, this would mean that the expressions themselves are equal. Also, it should be 
accounted that the function [p (z)] is an even one, and consequently only the even (non- 
positive) powers of z in the Loran decomposition of the two expressions should be taken 
into account. After performing the Loran decomposition, it becomes evident that equality 
of the two expressions is possible only if a = 4, b = 0, c = — 6OG4, d = — 140G6- Since 
these coefficients give exactly the algebraic equation (33), it follows that the Weierstrass 
function and its derivative (34) satisfy equation (33), thus representing "uniformization 
variables" for the equation (33), i.e. variables, which are functions of the complex variable 
z and at the same time satisfy the given algebraic relation. 
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It is important to stress that the "tilda" differentials dx 4 and dx 5 , which are related 
through the algebraic relation (33) and the parametrization (34) with the Weierstrass 
function, do not result in any dependence between the original differentials dx 4 and dx 5 , 
which should remain independent since are related to the independent coordinates in the 
gravitational Lagrangian. The reason for this independence between the tilda and the 
non-tilda differentials is that the linear-fractional transformation (31), which relates dx 5 
and dx 5 , introduces an additional arbitrariness in the non-tilda differentials due to the 
arbitrary complex functions a, b, c and d. 

In the present case, however, there are some specific facts about the parametrization 
of the curve: 1. The parametrization shall be carried out not in respect to two of the vari- 

~5 

ables, but in respect to dx and another variable, which is - - the two of the parameters, 
entering the linear-fractional transformation. The rest of the variables, entering the cubic 
algebraic equation shall be "hidden" in the free term, which is a function. So actually the 

- — '5 

final result will be for dx (or -), expressed through the Weierstrass function, but in order 
the parametrization to be consistent, the remaining four differentials (dx 1 , dx 2 , dx 3 , dx 4 ) 
should be related to a, b, c and d in a complicated way. 2. After performing the transfor- 
mation (31) with the purpose of chosing a, b, c and d in to eliminate the highest (third) 

- — '5 

power of dx ,the obtained equation will be like equation (33), but with g 2 and #3— func- 
tions and not complex numbers. On the other hand, the standard parametrization (34) 
with the Weiestrass function and its derivatives is valid only for g 2 and g^ complex num- 
bers. However, it will be proved in the next sections, that in such a case the formalism 
and the parametrizable equation can also be used. It will be shown that all the coefficient 
functions (those standing before the pole terms) in the Loran expansion can be found if 
the sums G n are known. 

V. TRANSFORMED CUBIC EQUATION 

WITH THE HELP OF THE 

LINEAR - FRACTIONAL TRANSFORMATION 

In order to derive this equation, all the terms with dx 5 in equation (20) shall be singled 
out and it can be written in the following way 

A{dx 5 ) 3 + B(dx 5 ) 2 + C{dx 5 ) + G (4 \dx 4 , ...dx^g^,^,^) = 0, (36) 

where A, B and C are the following functions, depending on g^, T^, Rij and the differ- 
entials dx a , dx 13 ; the indices a, j3 — 1, 2, 3, 4; r = 1, 2, ...5. 

A = 2pT r 55 g 5r (37) 

B = QpT r a5 g 5r dx a (38) 



and 

C = -2R a5 dx a + 2p(2T r a/3 g 5r + Y\ a g^ r )dx a dx fi . (39) 

The function G^(...) is of the following form 

G (4) (dx\ ...dx 1 , gij , Tj, Rik) = -R aP dx a dx p + pdx^ l dx a dx fi Yl t{a g^ r . (40) 

In (40) the indice 7 = 1,2,3,4 and (a,/3) means symmetrization in respect to the two 
indices. Further, after performing the linear-fractional transformation (31), one easily 
obtains the new cubic algebraic equation, written in terms of the new variables dx 5 : 

(G (4) c 3 + aQ)(dx b f + (bQ + aT + 3c 2 dG (4)) (^ 5 ) 2 + 

+ {aS + bT + 3cd 2 G (4) ){di 5 ) + (bS + G (4) d 3 ) = , (41) 

where Q,T, S denote the following expressions 

Q = Aa 2 + Cc 2 + Bac + 2cdC, (42) 

T = 2Aab + Bbc + Bad + 2cdC, (43) 

S = Ab 2 + Bbd + Cd 2 . (44) 

In fact, the linear- fractional transformation is performed with the purpose of setting up 

- — '5 

to zero the expression before (dx ) 3 , from where G^ is expressed as 

(?W = -^. (45) 

c 6 

This equation is the first additional equation, which is imposed in order to receive the 

parametrizable form of the cubic equation. Let us write down in more details equation 

(45) , in order to understand its meaning. Making use of the expressions for G*- 4 -* and 

Q, it can be written in the form again of a cubic algebraic equation in respect to the 

remaining four differentials 

pT r l{a g (3)r dx~ i dx a dxP + K { ^dx a d X P + K^dx a + 2 V (-) ' F r 55 g 5r = , (46) 

where K^l and K a are the corresponding quantities 

K% = -R aP + 2p- c (l + 2-) (2r> 5r + Vl a g Pr ) (47) 



KP = 2- 



and 

3p^T r a5 g 5r - (1 + 2-)i? a5 
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The indices a,/3, 7 = 1,2,3,4 ( but r = 1,2, ....,5) and (a, (3) means symmetrization in 
respect to the two indices. In other words, the imposed ("by hand") equation (45) simply 
fixes the cubic algebraic equation in respect to the remaining four differentials, 
if one would like to parametrize the differential of the fifth coordinate with 
the Weierstrass function. No ratios - and - are to be determined from this equation 
- later on from the equation in respect to the fifth coordinate they will be determined. 

~5 ~5 

Using expression (45), the functions standing before (dx ) , dx in (41) and also the 
free term in the same equation can be written in a form of an algebraic expression in 
respect to -, - and k 

bQ + aT + 3c 2 dG^ = d 3 {-3A-(-f + C-(-f + 2C~ - QC-+ 

c d d d dd d 

+3A 1& + B 111- 2B &- C 11^ (49) 

a d ddd d dd 

aS + bT + 3cd 2 G {A) = rf 3 {-3A(-) 2 ^ + B-(-f + 2C~- 

c d d d dd 

„aa „a „ab .a,b^ a, 

_ 3jB _ QC- + IB— + 3A- - 2 - C-}, 50 

cd c dd d d d 



and 



bS + GW = d 3 {-A(-f + A(-f + B(-f- 

c d d 

-B{-) 2 + C--C--2C-U. (51) 

c a c c c 

Let us now introduce the notations 

a ' — 5 

— = m dx = n . (52) 

c 

Equations (49-51) shall be rewritten in such a way so that the terms with powers of m 
will be singled out. The rest of the terms will be denoted by F, M and N and they will 
contain only powers of | and | only. The transformed equations (49-51), if substituted 
back into equation (41), allow one to write the equation into the following form 

-3A(-) 2 m 3 n 2 - 3A(-)m 3 n + [3A(-) 2 - r - 2B(-) 2 ]m 2 n 2 - 3B-m 2 n+ 
d d ddd d 

+ [-6C3 + B-(-) 2 - C(-) 2 ]mn 2 + [-6C + 2B-)- + 3A-(-) 2 - C^-]mn+ 
ddd d dd d d d 

+Fn 2 + Nn+[M- Am 3 - Bm 2 -Cm- 2-Cm) = . (53) 
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The terms F, M and N have the following form 

M = A( b 2 ) 3 + B( b 2 f + C- d , (55) 

In other words, we have transformed the original third-rank algebraic equation of five 
variables dx 1 , dx 2 , dx 3 , dx 4 , dx 5 into an algebraic equation of two variables only (m and 
n), but with a higher rank (in the case it's five). 

VI. A PROPOSAL FOR STANDARD 
PARAMETRIZATION OF THE CUBIC 
ALGEBRAIC EQUATION WITH 
THE WEIERSTRASS FUNCTION 

By standard parametrization it shall be meant that the cubic algebraic equation should 
be brought to its standard parametrizable form 

n 2 = 4m 3 — g 2 m — g 3 , (57) 

where g 2 and g$ are the already known complex numbers. Then one has the right to set 
up 

n = p (z) = — m — p{z) . (58) 

In order to obtain the parametrizable form (57), it is instructive to write down the obtained 
algebraic equation in the form of a third-rank polynomial of m with coefficient functions 
P\{n) , P2(n), Pz{n) and P±{n), representing quadratic forms of n and at the same time 
cubic algebraic expressions in respect to | and | : 

Pi(n)m 3 + P 2 (n)m 2 + P 3 (n)m + P 4 (n) = , (59) 

where 

P^n) = nn 2 + r 2 n + r 3 = -3A(^)V - ?>A-n - A , (60) 

P 2 (n) = qi n 2 + q 2 n + q 3 = [SAfy 2 ^- - 2Bfy 2 }n 2 - W-n - B , (61) 

a d d d 
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P 3 (n) = pin 2 + p 2 n + p 3 



d dd d 



n 



+ 



cb ,c ,b^ c 

-6C + 2B-- + 3A-(-f-C- 

dd d d d 



n-C-2-C 



P 4 (n) = Fn l + Nn + M . 
Let us write down the last expression in the following form 

1\ N9 M ( 1\ \"' 
2F J F \2F 

where n denotes 



~2 TT N 

n 2 + M - — , 
4F 



(62) 
(63) 

(64) 



n 



, . N 

F n + ^ 
2F 



In terms of n, the transformed equation (59) can be written as 

n 2 = Pi(n) m 3 + P 2 (^) m 2 + P 3 (n) m + P^n) , 
where the coefficient function P\{n) is 

Pi (n) = r i n 2 + r 2 n + r 3 = 



(65) 



(66) 



7-1-2 

■=n 
F 



and 



AT 
— s-ri - 


7-2 

l 


n + 


F 2 


i?2_ 





N 2 N 

'^1^=9 + r 2 ^ 

4F 2 2F 



r 3 



(67) 



- N — 

PAn) = ^-M 
y ' 4F 



(6£ 



The other coefficient functions P<i(p) and P 3 {n) can be written analogously, but with 
(<7i)92)<73) an d (pi,p 2 ,P3) in (67) instead of (ri,r 2 ,r 3 ). Note that unlike the expressions 
for r, q and p, representing cubic algebraic expressions in respect to j, and %, the corre- 



sponding "bar" quantities represent more complicated expressions, which are no longer 
polynomials. It is also not correct to consider the transformation from (p, q, r) to (p, q, f) 
as a linear affine transformation. The expressions iV and F, entering the coefficient func- 
tons of the transformation depend also on | and ^,so presumably they could also be 
expressed through (p,q,r). The above transformation shall be investigated further. 

Our purpose will be to identify the investigated equation (54) n 2 = Pi(n) 
m 3 + P 2 (n) m 2 + Pz(n) m + P^in) with equation (57) n 2 = 4m 3 — g 2 m — g 3 , for 
which we already know that the substitution (58) can be performed. In order 
to obtain the standard parametrizable form of the cubic equation, one has to require that 
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the two equations are to be made equal, which means that the polynomials Pi(n) , P^iji), 
Pz(n) and P±(n) (depending on the variable n) are to be made equal to the numerical 
coefficients 4, 0, —g% and —g 3 respectively. Therefore, the following system of equations 
should be fulfilled 

4 = Ti n 2 + r 2 n + r 3 , (69) 

= q x n 2 + q 2 n + q 3 , (70) 

-92 = Vi n 2 +p 2 n + p 3 , (71) 

AT 2 _ 
-* = W -M. (72) 

The last equation (72) represents the second additional equation, imposed in order to 
obtain the parametrizable form of the cubic equation. Note that this equation has an 
extremely complicated structure: since N, F and M are third-rank polynomials in 
respect to ^ and |, the equation will be of sixth order! This causes inconvenience 
in investigating such equations, therefore it is appropriate to search another variables, in 
terms of which the algebraic treatment will be comparatively more convenient. 

Let us try to find such variables. From the first and the second equations (70,71) the 
terms with n 2 can be excluded, and also from the second and the third equations. The 
obtained equations are 

4?i = (Mi - f{q 2 )n + (r 3 g x - riq 3 ) , (73) 

-mi = (Mi - PiQ 2 )n + (Mi - Ms) • ( 74 ) 

From the last two equations the terms with n can also be excluded and a fourth-rank 
algebraic equation is obtained in respect to pi, qi and r* (i — 1, 2, 3) 

(Mi - M2)( 4 ?i - Mi + riq 3 )+ 

+ (f 2 q 1 - f 1 q 2 )(g 2 q 1 + p 3 q x - p^) = . (75) 

The above equation represents the third additional equation, imposed in order to obtain 
the parametrizable form of the cubic equation. This equation is difficult to deal with, but 
there is a way to rewrite it in a more convenient and simple form. Let us introduce the 
"angular" type variables / and / with the corresponding components 

I — [l ,1 ,1 ) — (Z12, hs, hi) — 

= (M 2 - Mi, Ms - M2, Mi - Ms) , ( 76 ) 
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/—(/)/)/) — (/l2) /23, /3l) — 

= (ng 2 - Mi, ^2? 3 - f sg 2 > Mi - ng 3 ) . (77) 

In terms of these variables, the fourth-rank algebraic equation (75) will be reduced to the 
following quadratic equation 

Aq 1 l 1 +g 2 f 1 q 1 + l 1 f 3 + fH 3 = 0. (78) 

Having found the algebraic variety for (q 1 ,l 1 ,l 3 , f 1 , f 3 ), one can go back to find the al- 
gebraic variety for (p, q, r) . From there by means of the inverse transformation of (67) 

i 

T\ = —Fri ; r 2 = —F 2 r 2 — Nf\ (79) 





r 2 = —F 2 r 2 - 


- Nr 1 


N 2 


N _ 




4F n 


i f'2 f'd 

9P2 





r 3 = --=n =r r 2 - r 3 (80) 

Ah 2F 2 

(the same for (p,q)), one can obtain the "non-bar" variables (p,q,r). As already men- 
tioned, the coefficient functions of the above transformation depend in a complicated way 
on iV and F and therefore on 4 and %. Therefore, if the dependence of k and ^ on the 
"bar" variables is known, one would have a well-determined transformation (although it is 
a nonpolynomial and nonlinear one) from the "non-bar" variables to the "bar" variables 
(p,q,f). That is why the purpose in the next section will be to find this nonpolynomial 
transformation. 

Finally, it may be noted that if the parametrization of (57) with the Weierstrass 
function is performed and n — p (z) — -£■ and m = p(z), then the parametrized cubic 
equation can be written in an integral form 

dp(z) 



dz 



y / 4p 3 - g 2 p - g 3 p 

The variable z is a complex one and it may be different from the variables Z\ and z 2 , 
"performing" the complexification of the four- dimensional manifold (a = 1, ..4), and the 
integration is performed along some contour in the complex plane. 

VII. FINDING THE NONLINEAR 
AND NONPOLYNOMIAL 
INVERTIBLE TRANSFORMATION 

We shall start from expressions (60-62), from where one can find 



2 



.c r 2 



r 3 = -A r 2 = 3-r 3 n = -^ , (82) 

a Srs 
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3 ~d q3 
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9i = ~ r 2 ; + 2 

3g 3 d 9g 3 



J3) 



where it has been used that — 
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If expressions (83) for q = (q 1 , q 2 , g 3 ) are substituted 



into the defined by (67) expressions for q 1: it can be obtained 
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+ r 3 



3/V /c\ 2 6 



4^r 2 Vc?y d 



TV 



The corresponding equations for r = (ri,r2,r 3 ) are 

3r 3 /c\ 2 



r L 



c 1 
3g 3 -;---r 

W 172 



$4) 



36) 



37) 



TV 



r 2 



r 3 = n 



-3r 3 
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TV TV 
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F 2 
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From the first and the second two equations it can be obtained respectively 



N = —\-r 2 F 

ri 



fI 



and 



(3 + 4r 2 - 6f i 



-F 2 



(ri - l)r 2 



F+ .'2 



^3J 







39) 



(90) 



(91) 



12fi 2fi 4rf 

i 
where F = -^V- • It is important to note that Y can be found as a solution of the above 

d _ 

quadratic equation with coefficient functions, which consist only of r. Therefore 

1 



c 
d 



F* 



ZF 2 



N = OF i 



n 



-r 2 



1_, 



j_ 

F 2 



Now let us write down the corresponding equations for p from (62) 

d dd \d 



(92) 



(93) 
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(94) 



(95) 



If from the last expression C is expressed and is substituted into (93), an expression for 

P3ZF*[6 + ZF*] 



-, can be obtained in the form 

a 



b 
d 



Pi 



+ 



q 3 F 2 Z q 3 [2 + ZF' 



(96) 



where the derived expressions (92) have also been used. In order to obtain an expression 
for | m terms of the "bar" variables p = (pi,P2,P3), the "non-bar" variables pi,p 2 and p 3 
should be expressed from the system of equations for p : 



(97) 



Pi = 


_P1 

F 


N p 2 

; P2 = —Pi — r 

p 2 p 2 




P3 = 


N 2 N 
P \F* +P2 2F ^ 



(91 



and substituted into (96). The result is 



d q 3 



F 2 (6 + ZF 2 ) I 2 F 2 „ OF 2 Z -± 

— Pi ~ ~ ZJ 1 I — —ZPi + -^7— P2 + ZF*p 3 

(2 + ZF' 



(99) 



The only "unbar' variable q 3 can be expressed from the first two equations (84-85) for q 1 
and q 2 

* = -%* • < 100 > 

Also, from the third equation (86) for q 3 it can be obtained 

b 8 [QZq 3 + 30 2 Z 3 + (30 Z - 2) (Oq 1 + q 2 )\ 



3r.Q 



'd 30 2 Z 3 

and from the first equation! (3) the same expression can be found to be 



(101) 



3r 3 - 
d 



b 2Z 2 + q x 



Z 2 



(102) 



From the equality of the above two formulaes one relation between the 'bar' variables can 
be found. More concretely, since O and Z depend only on r, the relation will concern how 
g 3 can be expressed through O, Z and qi,q 2 - This will not be used further in the text, 
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since our main purpose will be to find the ratio l2 ^, which is to be used in the subsequent 
formulaes 



r 3 F _ F^(2Z 2 + q 1 )(2 + ZF l 
~qf' (Oft + g 2 )ATi 

where i^i is the expression 



(103) 



K l ^p 1 {2 + ZF^)-Z 2 {Q + ZF k ){ ( ^ + ^p 2 + p,) . (104) 

At this moment the only equation not yet used is the one, which can be derived from 
(93-95) for p 2 

2P3J(6 + §) , c r 3 1 r 3 3p 2 

r 3 3^(§) 3 (6 + ^ 

ft 2 (2 + §) 2 • (105) 

If the p- variables are expressed from (97-98) through their "bar' counterparts and all 

preceeding expressions are used, the following cubic algebraic equation in respect to 

i 

F 2 = T is obtained 

A^iT 3 + N 2 T 2 + N 3 T + iV 4 = , (106) 

where Ni,N 2 , N 3 and iV 4 are complicated expressions of the "bar" quantities only. These 
expressions will be presented in Appendix D. Therefore, the roots of this cubic equation in 
respect to T can be found and consequently, the quantities N, F, M, entering the second 
additional equation (72) also can be expressed in terms of the "bar" variables. This in fact 
proves that 1. The two additional equations (72) and (78), imposed in order to obtain 
the parametrizable form of the cubic equation, can be expressed in terms of the "bar' 
variables only. 2. The nonlinear and non-polynomial transformation from the (r,q,p) to 
the (r,q,p) variables is an invertible one. This is an important fact, since one first may 
study the properties of the algebraic equations, given by (72) and (78), and then chose 
the most convenient form for the ratios -, and §. 

a d 

VIII. PARAMETRIZATION OF A GENERAL 
CUBIC CURVE WITH COEFFICIENT 
FUNCTIONS OF A COMPLEX VARIABLE 

In this section an attempt will be made to deal with a cubic curve of a more general 
kind 

n 2 = M(z)m 3 + N(z)m 2 + P(z)m + E(z) , (107) 
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where M, N, P and E are functions of the complex variable z and therefore not complex 
numbers, as usually accepted in standard complex analyses 5 and algebraic geometry 13 . In 
other words, the main problem is whether it is possible to parametrize with the Weierstrass 
function the above equation, i.e. when does the Weierstrass function satisfy the equation 

A? = M ( z ) p 3 + N ( z)p 2 + p^ z)p + E ( z) ? (10g) 

dz 
As already briefly discussed in Sect. IV for the standard and usually investigated case of 
M, N, P, E - constants, the Weierstrass function parametrizes the cubic equation (107) 
only if M = 4, N = 0, P = — 6OG4 and Q = —140(76, but evidently in the present case 
of functions, the situation will be quite different. 

Let us first decompose p(z) into an infinite sum, assuming that | w | is a large number 
and therefore 

n _ j_ ^ r 1 i_ 

= \ + E A( 2 - + 3 (-) 2 + ••• + (" + 1)^ + -...) • (109) 

Z z z — ' ZD ZU ZU W n 

The first derivative of the Weierstrass function is 



dp 2 y-^ n(n + 1) 

p \Z) = 

and its square degree is 

4 ^ /M „ -■■; „_; V- " l"^U ,2(»-J.l ,||| 



P 



rf^ Z 3 ^^ Cc7 2 +™ V ' 



n(n + l) „, , v^n 2 (n+l) 2 



^2+n ^^ W 2 (™+ 2 ) 

ra=l n=l 



Note that in the strict mathematical sense, the second sum in the last expression is in 

00 00 

fact a double sum over m and n ^2 X] m " m +n — i-,2 m+n ~ 2 , obtained as a result of the 

ra=l n=l 

multiplication of the two infinite sums (110) for p (z) with different summation indices. 
Of course, since the two sums are equal and infinite ones, the representation in the form 
of a single sum is also correct. The appearence of the double sum should be kept in 
mind, since the idea further will be to compare the coefficient functions in the Loran 
power expansion of the functions on the left- and on the right - hand sides of (108), and 
naturally a double sum will appear in the R.H.S. of (108). 

In reference to this, an important remark follows. Suppose one works in the framework 
of standard arithmetical theory of elliptic functions, when M, N, P and E are assumed to 
be just number coefficients. Since the Weierstrass functions p(z) is an even function of 
the complex variable z (see Ref.5 and Ref. 7), the whole expression on the right-hand side 
(R.H.S.) of (108) will be an even one too. On the other hand, the function p (z) in the 
left-hand side (L.H.S.) of (108) is an odd one, but its square again gives an even function. 
Therefore, comparing the coefficients in front of the powers in z means that only the even 
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powers should be included in the infinite sum decomposition 



p'(z) 2 = 4-76G 6 -24G'A + ... , (112) 



z" z 2 

where G n will denote the following infinite sum of the complex pole numbers 

z — ' w n 

Note also another very important fact the proof of which is given in Ref.5: the infinite (in 
numbers of w) sum (113) is always convergent (i.e. finite) when n > 2, but for n < 2 the 
finiteness is not guaranteed! In the presently investigated case of M, N,P,E- complex 
functions, no information is available whether the R.H.S. of (108) is an even or an odd 
function in z. Consequently, one should not use formulae (112), but just start with the 
more general expression (111) for [p (z)] . 

In order to find the Loran decomposition of the functions on the right-hand side of 
(108), one should first find the second and the third powers of p(z), which may be written 
as 

i °° n-2 °° 2n 

P 2 W = 7Z + 2 E^+ 1 )7^r + 5> + 1 ) 2 ^ > ( 114 ) 

z *■ — ' W *■ — ' w 

n=l n=\ 

00 n-4 °° ^2n-2 



1 Z n ~ Z 

p 3 W = ti + 2 E^ + 1 )t^ + E^ + ia2 '' 



z* i —'~ ' W" '-^ ' w 2n 

n=l n=l 

00 n-2 °° 7 2n-2 °° Jin 

+E(-+ 1 )^r + 2 E(-+ 1 )^r + E(-+ 1 )V- ( 115 ) 

n=l n=l n=l 

Since these two expressions are to be multiplied by another infinite sums, here in (114-115) 
we have retained the single-sum representaion. 

The function E(z) has the following Loran expansion around the zero point 

CO CO oo , (0) 

m = J2 ^ m = E «^ m + E -^ > ( 116 ) 

m=— co m=0 m=l 

where oL and bm can be represented as integrals along some contour in the complex 
plane (w is a complex integration variable) 

og) = — / -^- 6g) = — f E(w)w m - l dw . (117) 

2m J w m - l dw m 2ni J K J v ' 

The coefficient functions in the Loran expansion of the functions N(z),P(z) and Q(z) 
will be denoted respectively by c m , cL and c m . Each term of the expression for the 
right-hand side of (108) is a product of two infinite sums, and the final result is 
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M(z)p 3 + N(z)p 2 + P(z)p + E(z) = JT Y,{c { S +4 + 2{n+l)G 



oo oo 

(3) 

m=— oo n=l 



+ {n + l) 2 G 2n C^ 2 _ 2n + (n + l)G n c2+ 2 -n + 2 ( n + 1 ) 2 G2nC ( l+2-2n + 

+ (n + l) 3 G 3 J*l 3n + c% +4 + 2(n + l)G n c%_ n+2 + 

+ (n+ l) 2 G 2n c^_ 2n + c« 2 + c£L n G n + c^}z m . (118) 

In principle, the general case for an arbitrary m may also be considered. Then the above 
expression should be put equal to formulae (111) for [p (z)] , where in the first sum 
one should set up 2(n — 1) = m and in the second sum n — 4 = m. In the first sum in 

(111) the summation will be over values of m = 0, 2, 4, 6, ...2k, and in the second sum 

over m = —3, —2, —1,0, 1,2, In such a case and for a given n, one would have to 

consider a recurrent (in n) set of seven algebraic equations in respect to the four Loran 
expansion coefficients c n ,oh, ,Cn ,oh, and for the seven values of m — —6, —3, —2, — 1, 0, 
2k, 2k + 1 (k = 1,2...). Therefore, the system of equations is predetermined, which 
enables one to find not only the the unknown variables (the coefficient functions), but 
also certain relations about the "coefficient" expressions, represented in the case by the 
sums G n . This is an important moment, which shall be worked out further in this paper, 
and indeed certain interesting relations will be found. Moreover, since the summation 
over m in the R.H.S. of (118) ranges from — oo to +oo, terms with values of m, different 
from the above written shall be present also in the R.H.S. of (81), but not in the L.H.S. of 
(111). Therefore, two additional algebraic equations may be obtained by putting m = —2k 
(k ^ 0, 1,3) and then m = -(2k + l) (k ^ 0, 1) in the R.H.S. of (118) and then setting up 
the whole expression equal to zero. In fact, effectively instead of two additional equations 
one may have just one additional equation by putting m = —k (k > 3 and k ^ 6), 
so the total number of equations will be eight. This complicated calculation for the 
general case has not been performed in the present paper, because due to considerable 
technical difficulties it would be impossible to reconstruct analytically the whole set of 
Loran coefficients c n ,c n ,c n ,c n as solutions of the above system of eight algebraic 
equations. However, the calculation will be performed in Appendix A for the simplified 
case, which will be described also below. 

In this Section we shall restrict ourselves to the case of negative-power expansion terms 
in the decomposition of [p {z)\ , obtained for values of m = —6, —2, 0, and the main mo- 
tivation for this is the analogy with the standard parametrization of the cubic curve. In 
the next sections the case of positive-power expansion will be considered too. Unfortu- 
nately, even under this additional assumption it is impossible to resolve analytically the 
corresponding system of algebraic equations, if some other simplifying assumption is not 
added. This assumption will be given in the next Section. 

The first recurrent relation for m = —6 is 
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4 = cS + 2(n + l)G n c^ ] n _ 2 + (n + l) 2 G 2n cS_ 2 n + in + l)G n cf_ 3 l_ n + 
+2(n + l) 2 G 2 JX2n + (n + l) 3 G 3n c m 6 _ 3n + cS + 2(ra + l)G n c£ 2 L4+ 

+(n + lfG 2 J%_ 2n + cL 1 ! + c { l\_ n G n + c^ . (119) 

For m = —2 the relation is 

-76G 6 = cf + 2(n + l)G n c^„ + (n + l) 2 G 2n cg 2n + (n + l)G n 41+ 
+2(n + l) 2 G 2n c£ 2n + (n + l) 3 G 3 „cLl + <S + 2(n + l)G„c (2 i +2 + 

+(n + l) 2 G 2 „cLl + 4 1} + cSG n + 4 0) . (120) 

The last relation for m = is 

-24G 4 = cf + 2(n + l)G„4-„ + (n + l) 2 C 2n c2 n + (n + 1)^^+ 
+2(n + l) 2 G 2n c2„ + (" + l) 3 ^^.^ + c 2 2) + 2(n + l)G n c (2 i+ 

+(n + l) 2 G 2 „cLl 2n + 4 1} + cS_ n G n + c { % . (121) 

To avoid the possible confusion why n appears in the R.H.S of (119-121) but not in the 
L. H.S., let us remind that the left-hand sides for [p (z)] in these three equations have 
been obtained by fixing both summation indices n (n = m) and also m (m — —6, —2, 0), 
while in the right - hand sides only the indice m is fixed and the indice n is left unfixed! 
What will be performed in the next Section will be for each value of m — —6, —2, to 
fix in an appropriate way the possible values of n. Therefore more than three algebraic 
equations will be obtained, in which there will be no summation left. 

From the above system of three recurrent algebraic equations (119-121), the infinite 
sequence of coefficient functions c„ \cn\cn and c n should be found and moreover, it 
should be proved that this sequence is convergent in the limit n — ► ±oo. Still, because 
of the restriction to three values of m only, even if is possible to find c n ,c n , c„ c„ , it 
would not be correct to assert that the Weierstrass function parametrizes an arbitrary 
cubic curve with coefficient functions of a complex variable . This problem probably may 
be resolved by means of computer simulations only. 

In the next section the system of equations (119-121) shall be used for parametrizing 
a more simplified cubic curve (without the quadratic in p(z) term). 
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IX. PARAMETRIZATION WITH THE 
WEIERSTRASS FUNCTION 
OF THE CUBIC CURVE 



p{z 



= V - 92{z)p - g3[z 



The form of the cubic curve is the same as the parametrizable cubic curve in standard 
algebraic geometry, 5,10,11,13,24 but here it will be with g 2 and g% - functions of a complex 
variable. The key problem, which can be raised is: does there exist an algorithm for 
finding out the sequence of coefficient functions in the Loran decomposition of g 2 (z) an d 
gz{z), satisfying the above algebraic equation, provided that its more simple form will result 
in the following restrictions on the coefficient functions of the already considered general 
cubic equation (108) 
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r (0) 7 m 



m=— oo 



(123) 



From the first sequence of equations one obtains for the coefficients c m and c 



-,(2) 



4 3) = 4 



£) = for all 



/;/ 



(124) 



^ =0 for all m ^ . 



„, • • , < ■ (125) 

Taking the above relations into consideration, equation (120) for m = —2 can be written 

as 



.(3) 



-24G 4 = 2(n + l)G n <%L n - %> _ Gn c%_ n - c^ . 



(i) 



-,(i) 



,(o) 



(126) 



For values of n = 2 and n — 1 from the above equation the following equations are 
obtained 

' (127) 



-24G 4 = 24G 2 - c { o ] - c { 1] - c W 4 G 2 - c m 2 , 



.(i) 



(i) 



,(o) 



24G 4 = c ( L> + G lC K _i + c^ • 
From the above two equations c_\ and c_l can be found 
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24G 4 + 24G 2 - c, 
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(129) 
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c (1) + c (0) 

c + c -2 



24G, 



From (89), the general recurrent relation for n = p > 2 can be obtained 
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24G 4 - c, 
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(130) 



131^ 



It is clear that for the determination of c_ A , c_$ and c[_ 2 _ p one has to know c_ 2 an d Cq . 
There is, however, one exception - in (129) G 2 = X] ~^2 may be a divergent sum, so then 
one has c_ A = 24 (since G 2 is in the denominator, when G 2 — » 00, the corresponding part 
of the expression will tend to zero). 

Further, for m = and keeping in mind (124-125), equation (121) will give 
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-76G 6 = 2(n + l)G„c^ n + (n + irG 2 „4-2„ + (n + l)G n c^l n + 



(3) 



.(3) 



.(3) 



.(1) 



JO) 



+2(n + l) 2 G 2 „c^ 2n - $> - c™G n - <$> . (132) 

(3) 
For values of n — 4, 2, 1, when there are non- vanishing values among the coefficients % , 

the corresponding equations are 



-76G fi = 40G, 
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-76G fi = 12Go - c 
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(134) 



-76G 6 = 48G 2 - 4 1} - c^JGi - c, 



jo) 



(135) 



The above linear algebraic equations can be solved trivially linear to find the coefficients 
c_4, c_ 2 and c 2 , which depend on c_[ and q, 
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(137) 



'138) 



c 2 ' - 76G 6 + 48G 2 - 4 u; - cLlCi . 

Note that these coefficients can be divergent if G 2 and GJ are divergent. Taking into 
account equation (138) and also (125) for the case m — 0, but for a general value of 
n = p 7^ 1,2,4, an expression for c_^, can easily be found 



Ji) 



,(i)/ 



fc _ (_48G 2 + c^Gx) 



(139) 
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This formulae should be compared to the previously derived formulae (131), setting up 
—2 — p = —k. From the two expressions c_ 2 can be expressed 

c<?l = 24G 4 - c« - %H(-48G 2 + c^G,) . (140) 

However, c_ 2 can be expressed also from the two formulaes (129) and (136) for c_ 4 : 



9i 

Gr4 



£°> = 24G 4 + 24G 2 - 4 1} + -^(48G 2 - c^G x - 40G 4 ) . (141) 



Comparing (140-141), an expression for c_[ can be found, which does not depend on any 
Loran coefficient functions 



(i) _ I6G2 [3G2Gk — G^Gk — 3G^Gk-2] 



— 1 /m l/m /nr /nr /nr 1 ' V / 

Substituting this expression into the formulae (139) for c_ k , one obtains the convergent 
expression (k > 2, k ^ 4) 



The obtained expression (142) for c_{ can be substituted into (140) to find a formulae for 
c_2, which will depend on Gk and only on the Loran coefificient function c 

c <«» = _ c « + 24G4 + 1 l^g_. (144) 

C?4 G2 

The above expression is well defined also when G2 — ► 00. It shall be proved subsequently 
that such a case will turn out to be impossible. 

Further, from (144) and expression (130) for c_ 3 it follows 



c 



(1) _ G4 lQGk-2 

Gk — TT-Gk-2 



•3 



G, 



(145) 



But since k in expression (143) can take a value k — 3, it follows also 

The comparison of the two expressions gives the following formulae for the infinite sum 

Gk '■ 2 _j 

G fc = 7G fe _2 = 7 s C fe _ 2s = = T^Ci for fc = 2p 

= 7 p Gi for fc = 2p + 1 , (147) 
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where 

7 = 2^-^. (148) 

Another recurrent relation for Gk can be found also from (145) 

n (Gfc_ 2 ) 2 (Gk-4) 3 (G k ^2s) S+1 / 1/|n x 

G k = —^ = 7^ ^ = ••• = 7^ • (149) 

This formulae for values oi k = 2p and fc = 2p + 1 , combined with the previous formulae 
(147), allows one to find an expression for G\ : 

2p — 1 p 

/"(2p— 3 fp—1 

G 1 = ^V = ^V • (150) 

/~<2p— 3 /^p-1 

( - J 2p Lj 2p+1 

T/ie Zas£ formulae is interesting, because it shows that the divergent in the general case 
quantity G± in the present case is expressed through convergent quantities only - G3 and 
G2 P +i{p is of course a finite number!). Substituting (150) into (149), one can get expres- 
sions for G 2p and G2 P +i : 

G2 P = G 3 I — j ; G 2p +i = G3 ( tt j • (151) 

It is seen that G2 is also expressed through convergent quantities. 
From equation (130) for m = —6 one obtains 

4 = cL 1 ] + cl 1 L n G n + c2 (152) 

(n = 1, 00). Since c_ 4 and c_g_ n can be found, c_ & can also be determined. It is clear 

that among the coefficients Cm two of them - c_ e and c_ 2 can be determined from (143). 
The other coefficients will be determined in the Appendixes. 

Let us summarize the obtained results in this Section and in Appendixes A, B and C 
by formulating the following 

Proposition 2 Let g2(z) and g${z) are functions of a complex variable, which have a 
Loran function decomposition g2(z) = Yl Cm z m and gz{z) = Yl Cm z m and satisfy 

m=— 00 m=— 00 

the algebraic equation [p (z)] = 4p 3 — g2(z)p — #3(2), where p(z) is the Weierstrass 
(elliptic) function. Then the following statements represent (only) necessary conditions 
for the fulfillment of the above equation: 

1. The poles of the Weierstrass function (even if they are infinite in number), must be 
situated in such a way so that the sums G\ = Y ~ ana " G2 = Y ~^ are convergent (i.e. 
finite). The sum G\ can be expressed through formulae (150). 
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2. All the coefficients c m and c m in the Loran positive- and negative- power expansion 
can be expressed uniquely from the finite sums G n . 

3. The sum G\ is proportional to the sum G 3 with a coefficient of proportionality, 
equal to the ratio of the sums G 2 and G 4; i.e G\ = qG 2 (from A22). This formulae 
follows also from the more general one G 2p +i = G^^r)^ 1 (151) for p = 0. 

4. As a consequence from the above relation and formulaes (147-148), the sum G2 can 
be uniquely expressed as G 2 = VG1G3 . 

5. All the even-number sums Gq,G 8 ,G w .... equal to zero. 

6. The following relation is fulfilled G\ = 20 2 k G 2 k+2 , which can be obtained from 

(B15) and (B16). In order this relation to comply with statement 5, additionally one 
should have that G 5 , G 7 , G13, G15, G 2 ±, G 23 .. should be zero. However, G 9 , Gn,Gn,Gi 9 
are different from zero. 

Finally, with the help of (152), a check can also be made for the consistency of the 
obtained results. Substracting the two equations (152) for values of n and n + I, one 
obtains 

a " Gn+l - = (153) 



(it's more appropriate to divide everywhere by G n ). From (151) for values of n = 1p and 
I = 2q, for example, it can be found 



G n +i 1 G 



M 



2q 



xr , ■ (154) 

For other combinations (even and odd) of n and I the calculation is similar. Using the 
above formulae, it can easily be verified that equation (153) is identically satisfied. This 
confirms that the investigated in this paragraph system of equations gives consistent and 
noncontradictory results. 

X. POSITIVE - POWER TERMS IN THE 
INFINITE SUM DECOMPOSITION OF 
EQUATION [V (*)] 2 = 4p 3 - g 2 {z)p - g 3 (z)- 
-THE CASE OF POLES NOT AT INFINITY 

For the purpose, the expansion (111) for [p (z)] has to be used, and a change in the 
summation indices is performed so that positive-power terms starting from n = 1 (without 
the free term) are taken into account 



„ CO 00 

p{z 



\~^/ -n2' ~n9^v 9n . X * 



n + l)\n + 2) 2 G 2{n+3) z 2n - 4 > (n + 4)(n + 5)G n+6 z n . (155) 



n=l ra=l 
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00 (1) 
Using the Loran decomposition g2(z) = J2 c m z m , it can be obtained for —g2(z)p(z) 

m=—oo 

for the positive terms only 

oo oo 

-g 2 {z)p{z) = -J2 <&*" -EE( n ^+ l)G n . k+2 c^z n . (156) 

n=l n=l k=—oo 

oo 

By means of the relevent expressions for p 3 (z) from (115) and for g 3 {z) = Yl Cm z m , 
the following expression can be obtained for the positive - power terms of 

2 



m=— oo 







p(z) -4p 3 + g 2 (z)p + g 3 (z) 



= ^2[A (n,G)z n + A 1 (n,G)z 2n + A 2 (n,G)z 3n ] , (157) 

n=\ 

where A (n, G), Ai(n, G) and A 2 (n, G) are the following expressions 

A (n, G) = 8(n + 5)G„+ 4 + 4(n + 3)G n+2 + 4(ra + 4)(n + 5)G n+6 - 



_ r W _ r (o) 

°n+2 c n 



^( n + A;H-l)G n+fe+2 c2, (15* 



fc=0 



Ax(n, G) = 4(n + 2) 2 G 2(n+1) + 8(n + 2)G 2(n+1) - (n + lf{n + 2) 2 G 2(n+3) , (159) 

A 2 (n,G)=4(n + l) 3 G 3n . (160) 

Two important observations can be made from these expressions. 

First, the unknown coefficient functions c n \ 2 , c n and c_ k are singled out only in 
A (n, G). Therefore, it is more appropriate to write down (157) in the form 

oo oo oo 

J2 M™, G)z n = -J2 Mn, G)z 2n - ^ A 2 (n, G)z 3n . (161) 

n=l n=l n=l 

The above mentioned coefficient functions can be determined in such a way so that the 
characteristics of the infinite sum on the left-hand side (L.H.S) would correspond to the 
characteristics of the infinite sum on R.H.S. Such a characteristic would be for example 
the " convergency radius" , defined in standard complex analyses (for the infinite sum on 
the L.H.S.) as 

R = lirrin^^ \ A (n, G |~« . (162) 
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Second, since the infinite sums in (158-160) are finite, it can be seen that 



lirrir. 



Ai(n,G 



UlTLr. 



A 2 (n,G 



oo . 



(163) 



Let us find the convergency radius of the first infinite sum on the R.H.S. of (158) 
Ri = liirin^oo | Ai(n,G \~™ = lim r 



In I Ai(n, G) I . 
,exp{ } 



n 



exp{—lim ri 



d\Ai\ 

dn 



■} 



(the Lopital's rule has been used, since n — »■ oo and | A\(n,G) \—. 

that 

d | Ai(n,G) | t 

2)G2(n+l) + 8G2(n+l) + 



(164) 
oo. It can be found 



dn 



[n 



+8G 2(n+1) - 2(n + l)(n + 2)(2n + 3)G 2(n+3) . 



(165) 



This expression is a third-rank polynomial in n, while | A\ \ is a second-rank polynomial. 

d\A,[ 

oo and consequently 



Therefore the ratio -rjV wn ^ evidently tend to infinity when 

Ri = exp[— oo] = . 



/; 



(166) 



The same can also be proved in an analogous way for R 2 . 

Now, since on the R.H.S. of (161) one has an infinite sum with a zero convergency 
radius, it's natural to suppose that the same holds also for the L.H.S. A simple calculation 
shows that this may happen only if 



limr 



lim r 



dc (2) 

"Si+2 I dc. 



d | A (n,G) 
dn 



(■» 



i 



dn 



dn &G n 



+ J2 Gn+k+2 C-k 



OO . 



rr 



-,(i) 



Jo) 



(i) 



This means that either c n l 2 , orcj, , or c_ k have to be proportional to 



(i+i) 

(i) 



(167) 
where / > 2. 



In such a way, we acquired information how the coefficient functions c n behave both in 
the positive-power (n > 0) and the negative power decomposition (n < 0). 



XI. POLES AT INFINITY (w -> oo) IN THE 



POSITIVE - POWER DECOMPOSITION 



OF 



p{z 



— /I ^3 



V - #2(»P - 9z(z 



If the period of the Weierstrass function can be represented asro = qw\ +p, then in the 
limit w — ► oo one has 



GL 



,G7 



^ ZJ7 n ^ foCCJi - 



(^l+p) T 






'168) 



where £(n) denotes the Riemann zeta-function. Therefore, in the asymptotic limit n — > oo 
(when this limit is used, as in the preceeding section), one should also take into account 
the asymptotic limit of the zeta-function as compared to the other power-like terms of n. 
For the present particular case, let us calculate for example R 2 - It can be found that 



lirrir 



and consequently, 



d\A 2 \ 
dn 



lirrir 



3 | C(3ra) - ra(ra + l)C(3ra + 1) 



n + 1 



C(3n) 



DC 



(169) 



In I A 2 I 
R 2 = lirn n _ 00 ex'p\ } = lim n ^ 00 ex , p\—oo) = . 



n 



;i70) 



After performing similar calculations, the same result can also be obtained for R\ . 
Let us find now Rq. Some more lenghty calculations will give 



lirrir 



d\A \ 
dn 



lirrir, 



C<U(n,k,p) 



, A I ,l ^\C( 2 \n,k,p) I ' 

where C^(n, k,p) and C^i^n, k,p) will be the following expressions 



(171) 



C {1) (n,k,p) 



d 



4(n + 4)(n + 5)(n + 6) l dn 



ii;lCn% + c r ° ) + J2an + k + 2)c^] + 



fc=0 



<(n + 7) + J2(n + k + l) 



fc=0 



dc {l) 
C(n + k + 2)—^ - (n + k + 2)c a l((n + k + 3) 

(III 



(172) 



C {2 \n,k,p) 



1 



4(n + 4)(n + 5)(n + 6) 



-,(i) 



.(0) , Y^/V, lJ.x1VLxI.xoU 1 ) 



C2 + d 0) + V(n + fc + l)C(n + fc + 2)c^ 



fc=0 



(173) 
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It is easily seen that Rq will again tend to zero if the nominator C^(n, k,p) in (172) tends 
to ±00. This may happen if at least one of the equations below is fulfilled 

dc^\ dcT 

— j-^— const.n ; " const.n (I > 3) , (174) 

an an 

, 00 
-f[J2^ n + k + 2 ) c -i] ~ constnr ( l > 3 ) 

=,"::./,,.., ._._,_.. ,„, 

00 00 1 

J2an+k + 3)c^ k =EE c -i^^ ~ constnr ( r>1 ) 

fc=0 fc=0 p^O 

==>- c_^, ~ const.n r p n+ +3 dn for some k and some p , (176) 

00 1 dc W 

^ ^ — Ti^—r 1 ~ const. n r (r > 2) 

^c a l~ const, f n r p n+k+2 dn . (177) 

It is interesting to note that after performing a change of variables 

tn = x exp(-t) = p , (178) 

the integral (175) for c_ k can be transformed as follows 
00 00 

' n l p n+k+2 dn= eXV[ ~ t t ^ 1 + 2)] I exp(-x)x l dx . (179) 

t 

Keeping in mind the generally known formulae for the Gamma function 
00 °° ™ 

exp(-x)x l ~ 1 dx = / exp(—x)x l ~ 1 dx + Y^ ui~^\ ' (^) 



n=0 



it can be derived for c_ k 



<£> " con S t.^m±^l{- / ^(-s^ds + T(/ + 1) - £ w , ( / f J • (181) 

"I n=0 ' 

Let us denote the integral in the above expression by 

t 
1(1) = / exp(-x)x l dx . (182) 



40 



Then the following relation can be found after performing multiple integration by parts 

1(1) = - | x l exp(-x) \ x Jz\ +1(1 - 1) = ....- | x l exp(-x) |*Z* - 

-/ | x^expi-x) \ x x =\ - l{l - l)..(l -k + 2)\ x l ' k+1 exp{-x) \ x x z\ + 

+1(1 -l)(l- 2). ..(I -k + 1)1(1 - k) . (183) 

Continuing in the same way, one can derive 

1(1) = -exp(-t)[t l + It 1 ' 1 + ..+ l\t) + 

t 

+eM - l)ll + l + lil - 1) + .. + P] + nJ eM - x)dx = 

1 

l-l , l-l 

= -exp(-t)^2— t l + exp(-l)^2l(l - l)..(l - k) - exp(-l)l\[exp(l - t) + 1] . (184) 

fc=0 fc=0 

Substituting this formulae into (181) and returning to the original variables (n,p), finally 
an expression is derived for c_ k : 

(!) ~ „™„w iM+i P n\f„ i „^f-\\\ , „,V ° +i 



c -k 



C0 ^-(- 1 ) i+1 7£vTI^ ! ^ + e ^- 1 )) + ^E^ 



(lnp) l+l1 '^ ' ry ,, v ^ dtk 

(-1)" 

n\(l + n + l) 



k=0 
l—l oo 



-exp(-i) £ /(/ -i)..(/- k) + r(i + 1) - E ^7rr^-TT > • ( 185 ) 



fc=0 n=0 

The above expression is of course not an exact one, since other representations of c k may 
exist, when the convergency radius of the infinite sum tends to oo. Nevertheless, it can 
be used as a possible model representation. 

It may happen also that the denominator C^(n,k,p) tends to zero, and then again 
the covergency radius will tend to infinity. For the purpose, the following two equations 
have to be fulfilled together 

ci+2 + ci 0) ~ const.rf (r = 1, 2, 3) (186) 

and 

C„ +fc+2 <££ - n =► c2 " const.n p n+k+2 . (187) 

The last possible choice is when the denominator is finite, but the nominator is infinite 
(equations 174-177). The denominator is finite when 

c n+2 ~ const.n 3 c_ k ~ const.n 3 (188) 

c W k ((n + k + 2)~ const.n 2 . (189) 

But since these equations contradict equations (175-176), such a case of finiteness of the 
denominator and tending to zero nominator has to be excluded from consideration. 
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XII. APPLICATION OF THE EQUATION 

2 

= 4p 3 -g 2 (z)p-g 3 (z)- 



p{z 
-ANOTHER WAY FOR PARAMETRIZATION 
OF THE CUBIC EQUATION (53) 

In Sect. VI a parametrization of the cubic equation (53) was proposed, based on presenting 
the equation in the form of a cubic equation in respect to one of the variables, and more 
importantly, applying the cubic equation for the Weierstrass function p(z) and setting 
U P 92, (73— complex numbers. In this Section it shall be demonstrated how one can obtain 
such an equation with #2, #3 - functions of a complex variable. 

It is instructive first to note that the last terms in the square parenthesis in (53) 
depend only on the variable n. Therefore, let us denote the term in the parenthesis by II 2 
and thus let us represent (53) in the form of two cubic equations 

n 2 = M- Am 3 -Bm 2 -C(l + 2-)fh (190) 

c 

and 

II 2 = P^njm 3 + P 2 (n)m 2 + P 3 (n)m + P 4 (n) . (191) 

The quadratic forms Pi(n)(i = 1..4) are the same as in (60-63), but with a reversed sign 
of n. Also, in the first three forms the last terms (—A), (-B) and (— C — 2-C) are absent. 
In order to bring the two equations to a parametrizable form, let us perform the linear 
transformation 

m = rm + s . (192) 

The first equation (190) transforms to the following equation 

h h 

IT 2 = [M - As 3 - Bs 2 -Cs- 2-Cs] + [-3rs 2 A - Cr - 2Brs - 2-Cr}m+ 

+ [-3r 2 sA - Br 2 ]rn 2 - Ar 3 m 3 . (193) 

In order to obtain the parametrizable form 

p(z)] 2 = 4p 3 -g 2 (z)p-g 3 (z) (194) 

of the cubic equation, one should require 

B , 4.i 2,4 x i 

S = ~JX r = (--£)* =i*(-j)*, (195) 

g 3 (z) = As 3 + Bs 2 + C(l + 2^)s -M=-M- ^C(l + 2^) , (196) 
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.2.2 , 2 



g J z ) = 3rs 2 A + Cr + 2Brs + 2-Cr = % -^-[C(\ + 2-) - — K . (197) 

a As d 3A 

Since in the preceeding sections it has been proved that the Weierstrass function p(z) 
parametrizes equation (194), one has the right to set up 

m = p{z) Il = p(z) (198) 

and consequently, the transformation (192) acquires the form 

m = rp{z) + s . (199) 

After performing this transformation, the second algebraic equation (191) assumes the 
following form 

P>)f = Qi{n)p\z) + Q 2 (n)p 2 (z) +Q 3 (n)p(z) +Q 4 (n) , (200) 

where 

Q^n) = P^ny , (201) 

Q 2 (n) = 3r 2 sP!(n) + r 2 P 2 (n) , (202) 

Q 3 (n) = 3rs 2 P^n) + rP 3 {n) + 2rsP 2 {n) = -g 2 , (203) 

Q 4 (n) = P 4 (n) + P 1 {n)s 3 + P 2 (n) s 2 + P 3 (n) s = -g 3 . (204) 

The first equation (194) and the second equation (200) will be identical in respect to the 
cubic and the quadratic in p terms if 

P 2 (n) B 4 4 , N 

s = i^- = r 3 = -^ = — - . 205 

3Pi(n) 3A Pl ( n ) A 

In fact, this means that the first two terms (cubic and quadratic) in the transformed equa- 
tions (201) and (202) are identical with the corresponding ones in the original equations 
(190-191): 

P^n) = -A P 2 (n) = -B , (206) 

Qi(n) =4 Q 2 {n) = . (207) 

In other words, the linear transformation is such that it transforms the straight line 

\h{n) + C(l + 2-)]m + P 4 (w) - M = , (208) 

c 
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defined on the points of the original cubic equation (53) (therefore - intersecting it), into 
the straight line 

(Q 3 (n) + 92(z))p(z) + (Q 4 (n) + g 3 (z)) = , (209) 

defined on the algebraic variety "points" of the transformed cubic curve (obtained from 
(201) and (206)) _ 

(Q 1 (n)-A)p 3 (z) + Q 2 (n)p 2 (z) + 

+(Q 3 (n) + 92{z))p{z) + (Q 4 (n) + g 3 (z)) = . (210) 

Note also that it cannot be assumed that 

Q 2 (n) = -92(z) Q A (n) = -g 3 (z) , (211) 

because then it will turn out that the linear transformation is a degenerate one and 
maps the straight line into the zero point. Of course, the choice of variables (206) does 
not mean that an additional and restrictive assumption has been imposed. The variable 
identification (206) simply helps to "fix" the transformation so that an one-dimensional 
submanifold (the straight line) is mapped again into an one-dimensional submanifold. In 
such a way, the original cubic equations (191-192) are replaced with the transformed ones 
(194) and (200), which differ yet in their last two terms. 

Now, using expressions (203-204) for g 2 (z) and g 3 (z) and also equations (205) for r 
and s, the following expressions are obtained, which will be used in the next section: 

b 1 4P 2 
2 (*)=4C(l + 2-W + -^-, (212) 



0=40(1 + 2^)1 


4/» 

+ 

3P? ' 


0(1 + 2^ 


- 2P| 
M I 



g ( z ) = -C(l + 2-)-^-M A- . (213) 

XIII. J - INVARIANT IN THE CASE OF 
ARBITRARY RATIO f AND 
IN THE GENERAL CASE 



The purpose in this section will be to see how the so called j or modular invariant of 
an elliptic curve, defined as: 

J(E) = 1728 -- 9 { 7 - 2 , (214) 



gi - 27s, 



3 



will change under some assumptions, for example when the ratio m — - of the parameter 
functions is an arbitrary one. In a broader sense, the idea is to see if there is any relation 
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between the possible motions and group transformations on the complex plane and the 
j— invariant. This will not be considered in this section. 

If the cubic equation (53) is satisfied under arbitrary rh, then all the coefficient func- 
tions in front of rh and its powers should equal to zero. Therefore, the following identifi- 
cation holds: 

M = h(n) -C(l + 2^)=P 3 (n) , (215) 

P^n) = -A P 2 (n) = -B . (216) 

Taking this into account, equations (196-197) for g 2 and g 3 can be written as 

2§ ' P 2 ~ 
<&(*) = — [-£■ ~ Ps] , (217) 

ps 3Pi 

g 3 (z) = P 3 5 - P 4 . (218) 

Pi 

Eliminating P 3 (n) from the two equations and expressing P±{n) gives 

92 P2 Po , N 

p 4 = £±_1 + _2_ _ ^ 219 

2§j5§ 3P 2 l ' 

Substituting P±{n) from (219) and P 3 (n) from (218) into equation (204) for Q 4 (n), and 
taking into account also (195), one obtains 

-g 3 = Q 4 ( n ) = -^ + S- + ^^-^2 • (220) 

27P 2 3 pa 

Performing the same substitutions in respect to equation (203) for Q 3 (n), a rather unex- 
pected result is obtained 

— 2§ P 2 2§ P 2 

-& = Q 3 (n) = — -£ -g 2 - —-^ = -g 2 . (221) 

6 p., 6 ps 

Therefore, g 2 for the cubic equation (190) does not change if assuming that rh is arbitrary. 
It follows also from (220) that g 3 = g 3 if and only if 

gl = -LJX . (222) 

y 2 g3 p4 V ) 

Since 

Pi = -A = -2pT; 5 g 5r and P 2 = -B = -Q p T r a5 g 5r dx a , (223) 
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clearly (222) will be a rather complicated algebraic equation of sixth rank in respect to 
the sub- algebraic variety of the variables dx a (a = 1, 2, 3, 4). 

However, it is more interesting to see when the j-invariants for the cubic equations in 
the case of arbitrary - — m are equal, i.e. 

](E)=j(E). (224) 

Making use of the defining equation (214) for the j— invariant and of equations (220-221), 
one easily finds 



8P| 2f P 2 
93 ~ Z7P? ~ ~^pl 92 



9t ■ (225) 



This is an algebraic equation on a Riemann surface. If one denotes 

g 3 (z) = w , (226) 

then in respect to w the surface has two sheaves - for +w and for —w. If the positive sign 
is taken, then the algebraic equation is satisfied for g 3 = g%— a case already considered 
(equations (223-224)). 

For the case of the minus sign —w, equation (225) will give another possible relation 
between the functions g2(z) and g${z) : 

APS 25 Po 
9s (z) = —L + T dr22 , (227) 

following from the equality of the j-invariants of the cubic equations for the case of 
arbitrary m. If the Loran decomposition of <?2(z) and g$(z) 



92{Z) 



E <& )zm 9*(z) = E c -^ m ( 228 ) 



is substituted into (227), then some additional relations may be obtained between the 
Loran function coefficients. 

One should keep in mind, however that the original cubic equation has been "splitted 
up" into two parts, so it should be seen how the j-invariants of the two parts are related 
to the j— invariant of the original equation. 

XIV. INTEGRAL REPRESENTATION IN 
THE LORAN FUNCTION DECOMPOSITION 
OF g 2 (z) AND A SUMMATION FORMULAE 
WITH THE FINITE SUMS G n 

In order to derive the summation formulae, let us remember that the coefficient functions 
Cq and Cm in the Loran infinite sum decomposition possess also an integral representation 
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c (l) = 



1 



92{W) 



2ni J tr 
c 



ra+l 



dw 



m = 0,l,2. 



(229) 



For the purpose of obtaining a formulae for g2(z), one should substitute the already found 
coefficient functions c_{ (eq.142), c__ 2 (eq.137), c_ 4 (eq.136) and c_ m (eq.143) into the 
Loran's decomposition formulae (228). The obtained expression is 
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(230) 



Let us see whether this will remain an expression for g2(z), or (which will turn out to be 
the case) the term g<i{z) on both sides of the equality will cancell out. For the purpose, 
let us rewrite the first two terms on the R.H.S. as 



ra=l m=l x 



qo(w) , \ „ 1 

^-rdw ) z m + 

m+1 / 2m 



w 



92{Wj 

w 



dw . 



(231) 



The integral and the sum in the first term can be interplaced, and also the formulae for 
the infinite geometric progression will be taken into account 



E 



1 



1 s 



So, expression (230) acquires the form 
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(232) 
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(233) 



according to Coushie's formulae. Therefore, g2(z) on the two sides of (230) cancels out. 
Changing the summation index in the last term in (230) from m to m = m — 3, one 
obtains the following expression for the last infinite sum in (230) 



00 , 



-,(1) 



m+2 



- 1 -z-' + 



cfi G, 



9 



TO=3 



I6G4 I6G4 G 2 4G 4 

3G 2 5 



)2T 2 + 



r< Ga rt 



^~ 3 + 



cL 1 ] Gl 



r<2 

Or 4 



20- I6G4 G. 



4 ^4 



(234) 
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where 



d m — _ g^ ■ (235) 



The essence of the above presented proof is the following: The Loran's decomposition 
(ranging from — oo to +00) is known to be convergent only in a segment. On the other 
hand, for the positive - power terms in to in Loran's sum one has the Coushie's formulae, 
where the integration is performed around a closed contour (a circle) in the complex 
plane. As for the other part with the negative-power in to terms, if z is not restricted in 
a segment area, in the general case there is no guarantee that the sum will be convergent. 
However, in the present case the negative-power part is represented by the infinite sum 
in (234), and also by terms with inverse powers of z, which are convergent in the limit 
z — > 00. The infinite sum J2d m z~ m is also convergent, since G m and G m+2 , entering d m , 
are finite when to — > 00 (according to a well-known theorem from complex analyses), and 
therefore the whole sum converges to zero in this limit (and also in the limit z — > 00). 
This is entirely consistent with the R.H.S. of (234), which is convergent to zero in the 
limit z — ► 00 and naturally does not contain a free term. 

XV. AN EXPRESSION FOR G\ FROM 
TAUBER'S THEOREM IN THE 
CASE OF INFINITE POLES 

Let us consider again the case of infinite poles, when G k — ► C(^)- Note that for k — 1 the 
zeta-function ((k) is not defined, so further the notation G± shall be preserved. 

Proposition 3 The coefficient function d m in the infinite sum (234) ^ n the limit to — > 00 
and under the assumption of infinite poles has the following behavior: 

dm = O(-) (236) 

m 

Proof. : ■ 

The proof is based on the following representation of the Riemann zeta-function, which 

can be found in Ref.40: 

11™ 1 

C(m) = + - + V B k m(m + l)....(m + k- 2) — - 

to — 1 2 L — ' k\ 

fc=i 



oc 



-m(m+l) (m + n-1) f BJx)x- m - n dx . (237) 

n\ J 
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Then d m can be represented as follows 



C(m + 2) - |Sc(m) 



iV _1 [l + 



C(4) 



% + l) N((2)(m-l) 



+ 



+ ^ ^7T7( m + 2)(m + 3)...(m + fc - 2) + ^^ (m + l)...(m + n - l)]" 1 , (23J 



fe=2 



where i^ , P m and iV are the following expressions 



Fk = (m + k — l)(m + k) — m(m + 1) 
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C(2) ' 



(239) 



— C(4) 

.B n (a;)a;~ m ~ n [-:r 2 (m + 2)(m + n)(m + n + 1) + -p-^-m}dx 



(240) 



N=l « 4 > 



2 2 C (2)' (241) 

In order to estimate d m , one needs the following inequalities for the last two terms in 
(238): 



B k F, 



J2 -4lj( m + 2)(m + 3)..(ro + fc - 2) > J] 



fc=2 



-Bfc Fk fc _ 3 



fc=2 



iV k\ 



^ p ip 

- — (m + l)...(m + n-l)>-— m 



(242) 
(243) 



From these inequalities, an inequality for d m also follows 



d m < N~ l 






EDk i'k 
~/V~iT 



(244) 



and (...) denotes the second term in the brackets in (238). It should be kept in mind also 

Ik 



that in the limit m — ► oo due to x m n 2 the ratio -^- is a very small term. 



Further, inequality (244) may be rewritten as 



d m < (^fm"- 1 )- 1 
n! 



m n m n l *-^ k\ Pm 



k=2 



(245) 



Since the last three terms in the squire brackets are small, the expression in the brackets 
can be decomposed (when X <C 1) according to the formulae 



[l+xr 1 = i-x+x 2 



(246) 
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Recall also from (240) that 



P„=0(±) F k ~rv?. (247) 

m 6 



Therefore, inequality (244) for d m can be rewritten as 

d m < n! (o(-Li) - 0(±) - 0(— 4 ) - J2 B^rn k - 2n A ■ (248) 

Neglecting all the small terms 0{.) and keeping in mind that in the last term in (248) 
the powers of m range from — 2n (k = 2) to — n — 2 (k — n), | d m |can be estimated 

| d m |< (n!) 2 J2 I ff I m ^ 2 ^ 2 < ( n! ) 2 ( n - X ) I ff ' m_n_2 = 

fc=2 

= 0(-4&) < (-) . (249) 

m n+z m 

This precludes the proof that d m = O(-). 

The purpose of the above proposition is to demostrate the opportunity to apply the 
Tauber's theorem in respect to the infinite sum (234), which for the presently investigated 
case and in terms of z — - can be rewritten as 

±Zr> = FM = {( a 9 ^< 3C(2) 



m=3 



16C(4)C(2) 4C(4) / V C 2 (4) 



- (1) ( c (1) 

' ^ + 77777777" + ^ 77771^-777 I "'} , (250) 



b c. 

'2C(4)~16C(4)C(2)^ ' V 16 C(4)G! '" Gl §| - C(3) 



where d m and c_i denote 

" m S C( m + 2)-SCN ' (251) 

_(i) = ^ , 48C(2)C(fc) ~ 16C(4)C(fc) - 48C(4)C(fc - 2) 
c -i"^ C (2)CW - C(4)C(* - 2) • (252) 

Let us remind the formulation of Tauber's theorem in Ref. 41 in terms of our notations: 

CO __ 

Theorem 4 // the infinite sum (250) ^2 dmZ 171 = F(z) converges to S when z — > 1 and 

m=3 

a/so d m = 0(— ), then the infinite sum J~J d m a/so converges to S. 

m=3 
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Applying this theorem to the infinite sum (250) and expressing Gi from there, one 
can easily obtain the following quadratic equation for G\ 



G\ C(2) 



16C(4)F 



~(i) 

16C(3)F + W) + 16 



Gl + ^J 2 P=0, (253) 

16FC 2 (4) 



where 



CO 



' = ^ 3 dm+ m-m- (254) 

Therefore, in the asymptotic limit z — > 1 and in the case of poles at infinity, Giis uniquely 
expressed through the Riemann zeta-function. As a solution of the quadratic equation, it 
has two values, which of course are finite, since the values of the zeta-function are finite. 
This again confirms the previously established result (although for a certain partial case) 
that Gi is a finite (convergent) quantity for the case of the investigated cubic equations. 

XVI. INFINITE POINT OF THE LINEAR - 
- FRACTIONAL TRANSFORMATION (31) 
AND THE TWO COUPLED ALGEBRAIC 
EQUATIONS FOR THE WEIERSTRASS 
FUNCTION 



In this Section again the obtained algebraic equation (53) will be studied, but for the case 
of the infinite point n = dx 5 = — - of the linear-fractional transformation (31), which for 
convenience shall be rewritten again 

~5 

, - adx + b . 

dx 5 = —^ (255) 

cdx + d 

. Substituting (255) into equation (53), a cubic algebraic equation is obtained for the 
variable m — - : 

c 

-Am 3 + 3A-m 2 + (B-, - A%- - 2-C)m + (-G (4) + 2-C + -C) = , (256) 

a d d 2 c c d 

where A, B, C and G 1 - 4 -* are given, as previously, by formulaes (37 - 40) in Sect.V. 

For the purpose of investigating the above equation, it is appropriate to remind briefly 
the two approaches, proposed in this paper. The first one was developed in Sect. VI and 
it contained in itself the following steps: 
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1. Transforming the original cubic equation into a cubic two- variable equation of the 
type (59), containing cubic terms of a new variable m — -, related to the linear- fractional 
transformation, and also quadratic terms of the original variable n = dx 5 . In a sense, the 
appearence of a new variable makes the equation more complicated, but at the same time 
this variable turns out to be the parametrizable one (with the Weierstrass function), and 
this shall be succesfully exploited also in this section. 

2. "Fixing up" the coefficient functions in the quadratic forms Pi(n), P2(n), P?,{n) 
(60-62), so that the standard parametrizable form 

n 2 = 4m 3 — g-im — g-$ (257) 

is obtained with g2,93~ definite complex numbers. The consistency of the parametriza- 
tion, in the sense of an additionally imposed condition, was ensured by the obtained 
quadratic algebraic equation (78) in terms of the "angular" type variables / = (li,h, h) = 
p x q and / = (/i, / 2 , / 3 ) = r x q 

±q 1 l 1 +g2f 1 <h + l 1 f + f 1 l 3 = 0. (258) 

Evidently, there is a certain indeterminacy and freedom in chosing the elements of the 
algebraic variety (I, /). However, by means of an appropriate choice of the ratio of the 
parameters in the linear-fractional transformation, equation (258) will acquire a very 
simple form, and in such a way this second difficulty can also be overcomed. 

The second approach, outlined in Sect. XII on the base of the analytical calculations 
in the preceeding sections, was based on parametrization with the Weierstrass functiion 
of the cubic form (257), but with g2 and g% - complex functions. In order to bring the 
original cubic equation into the form (257), a linear transformation with an appropriately 
chosen parameters was needed. 

In the present section, a combination of the two approaches shall be implemented. 
First, the original cubic equation shall be "splitted up" into the following two algebraic 
equations 

t 2 = -Am 3 -2-Cm+(-G {4) +2-C) , (259) 

c c 

o .b _ Q b 2 __ 6_ . 

t 2 = -3A-m 2 - A—m + B-m . (260) 

a d 2 d 

The first equation is a cubic one and the second - a quadratic one, and this fact is 

important. Moreover, there is no need to perform a linear transformation, and since it 

has been proven that an equation of the type (259) is parametrizable with the Weierstrass 

function, one can set up: 

t = p\z) m = p(z) . (261) 

Therefore, the second equation (260) can be written in the form 

p(zj\ 2 = Kp(z)[p(z)-a], (262) 
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where 

K = -?>A h - a= \l E. (263) 

d 3d 3A K J 

Introducing the notation 

g(z) = p{z) - a (264) 

and performing the integration in (262), one obtains 

dg* 



^Jl + f 



After two subsequent changes of variables 



Kdz . (265) 



i y 

g* =y and — = tgip , (266) 

one derives the following equation 

VKdz = 2 [ —dip = 2 [ -^— = In | ^- | , (267) 

J cosip J 1 — r 2 1 — r 

where 

r = simp . (268) 

In terms of the original variables, the solution is easily found to satisfy the equation 



2p — a + 2ep2 (p — a) 2 



a 



W/V**). (269) 



In (269) e = ±1 and Ai is the integration constant. Remembering the initial equation 
(262) and introducing for convenience the notations 



P(z) = ^y/Kexpi f VKdz) + ^- Q(z) = -UfK 



(270) 



the first order differential equation (268) may be rewritten in the simple form 

p\z)=P{z)+Q{z)p{z). (271) 

Performing again the integration, one obtains 

p(z) = ~ + H, (272) 

where H denotes 

H = -^exp( J Qdz) (273) 
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and A2 is the second integration constant. From (271) and (272) one can find 

^=7^rm and >'«=«*• < 274 > 

from where the integration constant A2 can be expressed, if one substitutes p (z) with its 
equivalent and well-known expression 

»'M = - 2 E77^W' < 275 > 



[Z — w 



In order to find the first integration constant, or at least a relation between the two 
integration constants, let us substitute p (z) from (272) and p(z) from (271) into the first 
cubic equation (259), which shall be written in a more general form 



p (z) 2 = B lP 3 (z) + B 2 p 2 {z) + B 3 p{z) + B 4 . (276) 



In the present case 



B 1 = -A B 2 = B 3 = -2-C B 4 = -G (4) + 2-C . (277) 

c c 

After the substitution, equation (276) acquires the following form in respect to the variable 
X — ^, containing the first integration constant Ai (which is to be found) 

B ± X 3 + B 2 X 2 + B 3 X + B 4 = , (278) 

where 

B 1 = -B 1 B 2 = ?>HB l + B 2 (279) 



B 3 = -3H 2 B 1 - 2HB 2 - B 3 B A = B 4 + B 3 H + B 2 H 2 + B X H 3 - Q 2 H 2 . (280) 

Note that the other integration constant A2 through H enters the coefficient functions 
Bi,I3 2 ,B 3 and 
B 4 . 



XVII. FINDING THE RELATION BETWEEN 
THE TWO INTEGRATION CONSTANTS 



The analyses in the preceeding section concerned the parametrization of the two coupled 
equations, and performing an integration only of the second one. Now we shall find 
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the relation between the two integration constants, inserting the found as a result of the 
integration solution into the first equation (276-277). 

Let us first observe that if one makes the formal identification 



A <& Bi 



B&Bs 



C <£> B ? 



D <£> G (4) , 



(281) 



then equation (281) is analogous to equation (36) from Sect. V. Therefore, one may apply 
the developed there approach in the same manner. 

After performing the linear-fractional transformation 



X 



aX + b 
cX + d 



(282) 



(where naturally a, b, c, d are different from a, b, c, d in the transformation (53 ,255)), and 
introducing the notation x — =, the transformed equation (278) is obtained in the form 



C 1 (X) X 3 + C 2 {X) X 2 + C 3 {X) X + C 4 (X) = 



(283) 



Since the subsequent step concerns the coefficient functions C\(X), C2{X), Cs(X), C^X), 
which are completely similar to the expressions (60-63) and (54-56), they shall be repre- 
sented below in terms of the new notations 



d(X) = r ± X 2 + r 2 X + r 3 = -3fii( = ) z X" - 3B 1= X - B 1 

d d 



(284) 



C 2 (X) = Ql X 2 + q 2 X + q 3 = [W^ift. - 2B 2 d) 2 ]X 2 - 3B 2 =X - B 2 , (285) 

a d d d 



C 3 (X)= Pl X 2 + P2 X+ P3 



— c — be — ,Cs 
-QB 3 = + B 2 = =-B 3 (=] 
d dd d 



X 2 + 



+ 



-QB 3 + 2B 2 l b = + 3B 1 ld) 2 -B T 



where in the last expression 



dd d d d 

C A {X) = FX 2 + NX + M 



X — B 3 — 2-B 3 
c 



(286) 

(287) 



— — b ,c „ — be 

F = B 3 = (=) 2 + 2B 3 = = 

d d dd 



(21 



c ,b , 



c b 



N = B 2 =(=Y + 2B 3 = = 



' d d 



dd 



(289) 
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and 

M = B~Af + B 2 {t? + b! . (290) 

add 

Now in turn is the folowing important assumption: In order to simplify the expressions 

(284-286) and to make the most appropriate choice so that a relation is imposed only on 

the coefficients a,b,c,d of the linear-fractional transformation (and not on B\,B 2 , B 3 , 

-B4, for example), let us suppose that C 2 (X),C 3 (X) are not quadratic, but linear in X, 

i.e. q 1= pi = 0. From (285) for C 2 (X) and (286) for C 3 {X) it follows 

b 2B 2 , , 

= = -=- , 291 

d 3B 1 K J 

c 2Bl , . 

= = -6 + _ A . 292 

d ZB l B 3 v y 

After performing the same change of variables as in (65), but in terms of our notations 

X = VF[X + ^=) => X = -^= - 4= , (293) 

V IF) ^F 2F ; 

the transformed equation (283) will be 

-2 



X = D 1 (X) X S + D 2 (X) X 2 + D 3 (X) X + D A (X) . (294) 

Let us find also the quadratic algebraic equation (258) for (Zi, l 2 , l 3 ) and (f\, f 2 , f 3 ) for the 
case pi — qi — 0. We have 

I = (h, h, h) = (0,M3. °) / = (/i» /2, /s) = (nq 2 ,r 2 q 3 , -f{q 3 ) . 

(295) 
It can easily be checked that eq.(258) is identically satisfied! In other words, we have 
succeeded to find such a (trivial!) algebraic variety (p, q, r), so that the quadratic algebraic 
equation (221) holds and therefore, the parametrization with the Weierstrass function can 
be applied: 

X = p'(v) x =^=p(v). (296) 

c 

In principle, v should be another complex variable, but here for simplicity it shall be 

assumed that v — z. Now it is interesting to note that relations (291) for =,(292) for = 

and (296) for = allow us to determine all the ratios between a, b, c, d. Consequently, F 

and iV in (288 - 289) are also determined, and taking into account transformation (293) 

and expresion (274) for A2 (expressed through p (z)), the linear- fractional transformation 

(282) can be represented as 

x = jJL+J = p(z)[F^X 2 exp(jQdz)-§} + l ^ 
x + i F~h 2 exp(jQdz) + l 
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As already mentioned, the ratios b : c,d : c can easily be found from (291-292) 

b _ B 2 B 3 d _ 3B 1 B 2 

c b\ - 9B{B 2 Z 2(B~l - W{B 2 ) 



(298) 



On the other hand, X = ^ can be determined also from (270) 

X = — = -eX^xpi / VKdz) - ea . (299) 

From the two equations (297) and (299), the first integration constant Ai can be expressed 
through the other constant A2. Moreover, (297) can be substituted into (292) and after 
some transformations, the following quadratic equation can be found for A2 

\ 2 2 + K 1 \ 2 + K 2 (p(z)-=) = 0, (300) 

d 

where _ 

K x = -^-=-exp{^- I VKdz)[2F -N^] , (301) 

2v F c 2 J d 

K 2 = -VfiVKexp(e f v^Kdz) , (302) 

L C J 

The two roots of the quadratic equation (300) can be expressed as 



K 1+ eJK 2 -A(p(z)-^K, 



<i> 



2 



A 2 = * a - . (303) 

But A2 can also be expressed from (274) and (262) in the following way 

Aj= VM]*_ = ^»-y. (304) 

exp[— e J yKdz\ exp(—e J yKdz) 

Now taking the square of A2 from (303) and expressing the equality of (303) and (304), 
one can obtain 

2Kp 2 (z) + 2(K 2 - aK)p(z) -L = 



where L denotes 



-eK x exp{-e I \fKdz)\lKl - A(p(z) - =)K 2 , (305) 

I = 2K 2 = + exp{-e I VKdz)K 2 . (306) 

d J 
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In order to get rid of the square root in (305) and thus to obtain an algebraic equation, let 
us take the square of both sides of (305). The result is a fourth-order algebraic equation 
in respect to the Weierstrass function p(z) 



where 

E x = AK 2 



E lP \z) + E 2 p\z) + E 3 p 2 (z) + E A p{z) + E 5 = , 



(307) 



E 2 = 8K(K 2 - aK) 



E 3 = 4[(K 2 - aK) 2 - LK , (30J 



E R 



E 4 = 4[K 2 K 2 exp{-2e / VKdz) - L{K 2 - aK)} 



L 2 + K*exp(-2e f VKdz) + 4=K 2 K 2 exp{-2e / VKdz) 



(309) 



(310) 



In a sense, it is surprising that (307) is a fourth-rank algebraic equation, while after 
parametrizing the "original" system of two coupled equations (259 - 260), one is left with 
a third-rank algebraic equation. It can be rewritten as 



P 3 W 



o b 2, 



, ,Bb b 2 Cd, . N 



G^ b C 



cA 



A 



dA' 



(311) 



It should be remembered, however that the fourth-rank algebraic equation (307) has been 
obtained after the differential equation (262) has been solved, and through the integration 
constant this amounts to a redefinition of the complex variable and the functions. Also, 
a more "specific" form of the linear-fractional transformation (282) has been performed, 
which in a way 'adjusts" the parametrizable variable in (294) X = p (v), and it was 
assumed that v — z. What is perhaps interesting, from a purely mathematical poini of 
view, is that this "redefinition" and transformation resulted in a higher-rank algebraic 
equation. 

It is important to note that in writing down the system of the two coupled equations 
(259-260), only the proven in the previous sections fact about the parametrization of a 
cubic equation of the kind t 2 = 4m 3 — g 2 {z)m — g${z) has been used. Since (t,m) turn 
out to be parametrized with (p (z),p(z)), and (t,m) enter also the second equation (260), 
this explains why the Weierstrass function "parametrizes' the second equation. But this 
does not mean yet that the more general cubic equation (256) can be parametrized with 
the Weierstrass function. As it has been already mentioned in Section VI, this fact needs 
to be proved, if this is possible at all. In the next Section, some necessary conditions 
will be found for parametrizing the cubic equation (256) of a more general kind (with 
all non-zero coefficient functions) . Although for a particular equation, the investigation 
in the next section might provide a partial answer to the raised in Sect. VIII problem 
whether it is possible to parametrize a cubic equation in its general form. 
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XVIII. FINDING THE NECESSARY 
CONDITIONS FOR PARAMETRIZATION 
OF THE CUBIC EQUATION (256) 

The analyses in this section concern equations (307) and (311). If p 3 (z) is substituted 
into equation (307), then the following equation is obtained in respect to p(z) 



Sip\z) + S 2 p 2 {z) + S 3 p(z) + S 4 = 0, (312) 



where 



(K 2 - aKf + 6-K(K 2 - aK) 



S x = AK 2 S 2 =A 

S 3 = 4{K 2 K 2 exp(-2e I Kdz) - L{K 2 - aK) + 



(313) 



s ,Bb b 2 Cd x , x 

+ 2K(K,- a K)(-- r - 2 -2--), (314) 



-K\ex V {-2e I \ficdz) - 4:=K 2 K 2 exp(-2e I \T~Kdz) 



(315) 



Two possibilities result from the above equation (312), and they shall be investigated 
here. 

The first possibility is when (312) holds for arbitrary values of z and p(z) , and therefore 
the coefficient functions S\, S 2 , S3, S4 must be identically zero. This means that 

5 1 = K = . (316) 
But if K — 0, then from the defining equality (302) for K 2 it follows 

K 2 = . (317) 

. If (316-317) are satisfied, then it will follow also 

5 2 = S 3 = , (318) 
and this is exactly what is needed. Now it remains only to impose the condition 

= S 4(K=0) = -I 2 - K\ = -1K\ . (319) 

Using expressions (301) for K\ and (288-289) for F and N , it can be shown that the above 
equation is equivalent to the following equation 

: . (320) 



b 




— cb 




AB 3 - 


-B 2 = = 


1 




dd 
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The first opportunity for this equation to be satisfied is when 

b 2B 2 2HB 1 . , 

= = -J^ = ^=± = . 321 

d W x 3Si 

Since we do not want to impose the rather restrictive condition B\ = —A = 0, it remains 
(with account also of (273)) that 

2A e r 

= H = -L e xp(-- / \T~Kdz) =► A 2 = . (322) 

s\jK 2 J 

If A2 = 0, then from (297) and (299) it can be found that 

Ax = -a + ep(z) . (323) 

But perhaps it is more interesting to see from (270)) and (271) that for the particular 
case it will follow 

P = Q = p'(z)=0 =>• p{z) = const. . (324) 

This represents the first found very nice and simple condition under which a parametriza- 
tion of (256) with the Weierstrass function is possible. The second possibility for equation 
(320) to hold is when 

4 _ _ _ c5 = iwiEu^m-^] = o 

Additionally, if equations (277), (279-280) and (273) for H are taken into account, the 
above equation can be represented in the following form in respect to the integration 
constant A2 

-16.3 3 A 2 A2 = . (326) 

Again, since A ^ 0, we come to the case A2 = 0. 

XIX. ANOTHER NECESSARY CONDITION 
FOR PARAMETRIZATION OF (256) 
AND APPLICATION OF THE THEORY 
OF RIEMANN SURFACES 

Now the second possibility when equation (312) holds shall be investigated. If one denotes 
p(z) = w, then the equation 

S lW 4 + S 2 w 2 + S 3 w + S 4 = (327) 
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and the couple of complex variables (w, z) determine a Riemann surface for the equation 

w = F(z, S^z), S 2 (z), S 3 (z), S 4 (z)) , (328) 

which represents a solution of the complex algebraic equation (327). Since the equation 
always has complex solutions and the Weierstrass function in principle is also a complex 
function, let us write w as 

p(z) = w = W\(z) + iw 2 (z) . (329) 

Remember also that no matter that from the original definition of the Weierstrass function 
in (35) it is not evident that the "imaginary" part tu 2 ( z ) is present, this can be guessed 
from the presence of complex poles, which after decomposing into real and imaginary 
numbers will give the "wi(z)" and "w 2 (z)" terms. If one substitutes (329) into (327), 
then equation (327) for the Riemann surface "splits up" into two equations for a couple 
of Riemann surfaces, determined by the complex pairs (wi, z) and (w 2 , z). The obtained 
equations for the real and imaginary parts of the complex algebraic equation (327) are 
respectively 

S^wf - w 2 2 f - ASxwlwl + S 2 {w\ - w 2 2 ) + S 3Wl + S 4 = , (330) 

w 2 [4S lWl (wl - w 2 2 ) + 2S 2Wl + S 3 ] = . (331) 

The last equation for the imaginary part is identically satisfied when w 2 = 0. Then the 
corresponding w± has to be determined from (330), which becomes a quartic equation 

S lW t + S 2 wl + S 3Wl + S A = (332) 

with roots x±, x 2 , x 3 , x^, satisfying the relation — x 3 — rr 4 = x\ + x 2 = u, and according to 
the general theory 13 the roots can be found after solving the equation 

- 6 + 2| W 4 +[(|) 2 -4K-(|) 2 = 0. (333) 

Dl 0\ 0\ Dl 

The other case when (331) is identically zero is when the expression in the square brackets 
is zero, and from there 

" 2 = 4S^ • (334) 

Again, substituting (334) into (330) for the real part and setting up w\ = Wi, one can 
obtain the following third-rank equation for W\ 

QASfwf + 32 1 S 1 2 ,S 2 wi + (4S1S 2 - 163fSi)wi - S^ = . (335) 

This equation has three complex roots for w\ and therefore - six roots for w\. Note also that 
each time when e appears in the coefficient functions, then for each function S±, S 2 , S3, S4 
and for every root, the corresponding expressions for the roots have to be taken once with 
e = l and once with e — — 1. So there will be a finite number of combinations, but more 
than six. 
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XX. NECESSARY AND SUFFICIENT 
CONDITIONS FOR PARAMETRIZATION 
WITH A CONSTANT WEIERSTRASS 
FUNCTION AND THE RESULTING 
NONLINEAR EQUATIONS 

In this section it shall be assumed that a constant Weierstrass function p(z) = const = 
w = W\(z) +iw2(z) parametrizes the algebraic equation (327) , defined over the Riemann 
surface (w,z). If the Weierstrass function is constant then it follows that 

dw dw 

-d~z = -d~z ' (336) 

and therefore 

dw! .dw 2 = dw! .dw 2 = , . 

dz dz dz dz 

Consequently, 

£-£"■ < 338 » 

Taking this into account and differentiating by z equation (334) for w 2 , one can obtain 

X' 

Wl = -yr, (339) 

where X and Y denote 

X -§1 y -f- < 340 » 

and ' means a derivative in respect to the complex variable z. Denoting f3 — w\ — const 
and inserting back w\ from (339) into (334), one can obtain the following nonlinear 
equation 

(X') 2 Y' 

\tY~ x x 7 + y ~ (3 = - (341) 

Introducing the notation Z = yr allows us to rewrite the equation as a cubic algebraic 
equation in respect to Z 

Z 3 + (Y - 13) Z - X = . (342) 

Finding the roots of this equation, one can easily express X as 

X' = F(X,Y)Y' , (343) 
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where Fi(X, Y) is the expression for the first root of the cubic equation 



The other roots are correspondingly 

Ux, Y) = -jEZM - £|S , ,345, 



Further the calculations will be done only for the first root Fi(X, Y). 

Now let us perform a differentiation by z of equation (330) for the real part, taking 
into account also the defining expressions (340), (339) and introducing also the notation 
T = Jr-. The result be another nonlinear differential equation 

2YY' - 4^—!- + T' = . (347) 

Equations (341) and (347) represent a system of coupled nonlinear differential equations 
for X and Y. Now we are going to establish an interesting mathematical property of this 
system : due to the peculiar structure of equation (341) as an algebraic equation with 
respect to the ratio of the derivatives of X and Y and subsequently expressing X as 
linearly proportional to Y (with a coefficient function F depending only on X and Y), 
this system will turn out to be an integrable one! So, we shall treat equations (343) and 
(347) instead of (341) and (347) . Substituting (343) into (347), we receive 

2YY' - AF 2 Y' + T' = . (348) 

Fortunately, this equation can be integrated to give 

T + f(X, ....) = -Y 2 + ^i§(Y - (3) 2 + Ah + 4/ 2 (349) 

and the integrals can be written as 



■■', 



h,h= I \^\/^+ {Z w l ) dY ■ ( 35 °) 

The integration function f(X, ...) depends on X and other constants or functions, different 
from Y. Now it remains to calculate (analytically) the remaining two integrals. For the 
purpose, let us perform the following three consequent variable changes 

^ = t , (351) 
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:^)*t (352) 



and 



^V / rT^= V /l + A^^. (353) 

It is straightforward to check that the integrals can be brought to the following form 

h = J{l + t)n<R=h{l+t)l \IZ\ ~J(l + t)ldi , (354) 

h = [(1 - t)Hdt = --t(l - t)f \ z = z z \ -- [(1 - t)i<M . (355) 

Note that if a closed contour is chosen, then the first terms in the above two integrals will 
be zero. If the integrals are calculated and expression (353) is used, then the following 
expression is obtained 




h + l2 = --rz< l + \ l + -^ v " 7 + 1-4/1 + —- ^ -Tl \. (356) 



Substituting back into expression (349), one obtains a solution of the nonlinear differential 
equation (348), which depends on terms in powers of fractional numbers. 



XXI. CONCLUSION 



Let us summarize the obtained results. 

In this paper a cubic algebraic equation has been obtained in respect to the differentials 
dX l of some generalized coordinates X 1 . The derivation of the equation was possible due 
to the representation of the contravariant metric tensor in terms of differential quantities. 
Also, in Sect. Ill the equation was derived upon assuming that dX l is either an exact 
differential, or that dX l are zero-helicity vector field components. 

The derived equation (20) clearly reflects the structure of the gravitational Lagrangian, 
and can be regarded as an equation for its all possible coordinate transformations (ad- 
missable parametrizations), provided the Christoffell's connection T^- and the Ricci tensor 
Rij are given. 

The main problems, which one encounters when investigating such algebraic equations 
are several, and in this paper only one of them is resolved in more details. 

The first and most serious problem is that the equation is defined on an algebraic 
variety of several variables, since in gravity theory one usualle deals with at least four- 
dimensional (and higher-dimensional also) manifolds. At the same time, the standard and 
known methods from algebraic geometry for parametrizing algebraic curves by means of 
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the Weierstrass function concern only algebraic curves of two variables. That is why in 
the paper one of the variables - dx 5 has been singled out ot the base of the physical 
consideration of the Randall- Sundrum models, and the other variable for convenience 
is chosen to be the ratio - of the functions a(z) and c(z), which enter in the linear- 
fractional transformation of dx 5 . The rest of the variables dx 1 , dx 2 , dx 3 , dx 4 enter the 
cubic equation in scalar quantities (functions). Of course, if a two-dimensional manifold 
is considered, then the variables might be related to the two variables of the manifold. 
Such an analyses of a two dimensional algebraic equation and its parametrization may 
find application in string and brane theory (also in gravity theory) . 

For the purpose of higher-dimensional algebraic varieties and equations, probably 
methods from the theory of abelian varieties and hyper-elliptical (Weierstrass) functions 
have to be applied. However, these methods are developed at an abstract mathematical 
level, 42 far from being adjusted to any concrete application. 

The second problem concerns the methods for bringing the algebraic equation to a 
parametrizable form. The standard approach of applying a linear-fractional transfor- 
mation has been chosen for the purpose, but in Sect. XII it was demonstrated also 
how a linear transformation can also be applied. In the last case, a couple of cubic 
equations is investigated, and what is interesting is that by means of parametrization of 
the first equation with the Weierstrass function p(z) the second equation is obtained in a 
parametrizable form, but in the general form of a cubic equation in respect to p(z). This 
justifies the performed investigation in Sect. VIII, concerning the (eventual) possibility 
for the Weierstrass function to satisfy a cubic equation of the general kind (108). Sect. 
XVI - XIX demonstrate how a differential equation with the Weierstrass function can be 
integrated, and after performing an additional linear- fractional transformation and find- 
ing the relation between the two integration constants in Sect. XVII, one comes to a 
purely algebraic (and not differential) equation. Interestingly, in some special cases, for 
example the case of parametrization with a constant in Weierstrass function in Sect. XX, 
one can effectively work with the approach of Riemann surfaces and even to obtain an 
integrable equation at the end. It should be noted here that this is possible in some cases 
only. For example, if one changes the formulation of the problem and would like to find 
the conditions for parametrization with a real-valued Weierstrass function, then it may 
be shown that the equation will not be an integrable one. This is not performed in this 
paper. 

As for the linear-fractional transformation, its advantages are the following: 1. 
It contains more parameters (in the case - functions a, b, c, d of a complex variable) and 
it makes possible to take account of the point at infinity. In principle, in complex 
analyses and projective geometry this is a well-developed procedure, but as far as physical 
applications in gravity theory and in relativistic hydrodynamics are concerned - this is 
still un unexplored area. 2. In Sect. VI it was proved that by means of a suitable change 
of variables it is possible to derive a second-order (quadratic) algebraic equation 
in terms of "angular"- type variables from the initial cubic algebraic equation. Since 
a quadratic equation is easier to deal with, this simplifies the analyses and moreover, a 
transition to the original variables can also be performed. 
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The third problem, investigated in more details in the present paper, is the form 
of the parametrizable cubic curve. This was discussed in the Introduction, and 
evidently a concrete physical problem from gravity theory has shown the necessity to 
investigate the case when g 2 and g% are complex functions and not complex numbers, 
as it is in the standard theory of elliptic functions. In regard to this, in Sect. VIII 
and Sect. IX two mutually related with each other problems for resolving are being 
stated: 1. Can the Weierstrass function p(z) parametrize an arbitrary cubic curve 
with coefficient functions of a complex variable ? 2. Can the Weierstrass function 
parametrize the well-known parametrizable form of the cubic equation, but again with 
coefficient functions depending on a complex variable? Although the explicit form of 
the equations for the Loran coefficient functions are presented in Sect. VIII, the first 
question still remains unanswered, and perhaps computer simulations only can help for 
its resolution. As for the second question, the answer is affirmative, and after solving a 
system of algebraic equations for various values of m and n, the explicit form of all the 
Loran coefficient functions c rn , Cm was found, both from the negative- power and the 
positive- power expansion. A confirmation of the consistency of the derived equations is 
equation (153) for a value of m = —6, which is being satisfied by the previously derived 
equations. However, the values of the coefficient functions are perhaps not so important 
as the result, which follows from this calculation, namely: the infinite sums Giand 
G 2 , which in the general case might be divergent, in the particular case of 
the "parametrizable" form of the cubic algebraic equation with g 2 = Qi{.z) 
and g 3 = g 3 {z), should be convergent! This fact, although of pure mathematical 
nature, probably deserves more attention and further elaboration from another point of 
view and by applying different mathematical approaches. The finiteness of G\ and G 2 is 
not imposed "by hand", but is obtained as a consequence of the fulfillment of the above 
mentioned equation with the Weierstrass function. 

It is to be noted also that if some other assumptions are made - for example in Sect. XI 
about poles at infinity in the positive-power decomposition of g 2 {z) and #3(2), then within 

the large n asymptotic approximation, a method for comparing the convergency 

radius of two infinite sums may be used for the determination of the coefficient functions 
c_ k . Comparing the expressions for c_ k in the negative- power Loran decomposition and 
in the other case of positive- power Loran decomposition with poles at infinity, it is seen 
that in the two cases c k are expressed in different ways. For example, in the first case 

from (143) cj k is expressed through the sums Gk, while in the second (asymptotic) case 

G n do not appear (185), and instead of them the Gamma function T(l) appears and some 

i—i 00 

finite and infinite summation formulaes of the kind ^ /(/ — 1 )...(/ — /c) and Yl n \(i+i+n) • 

fc=0 n=0 H ' ( n> 

Therefore, the assumption about poles at infinity and the consequent appearence of the 
zeta-function substantially changes the calculations. However, while the infinite point of 
the linear- fractional transformation may have some physical justification, the "poles at 
infinity" case for the moment does not have a definite physical meaning. 
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APPENDIX A: ADDITIONAL SYSTEM OF EQUA- 
TIONS FOR m= -3, -1 

In this Appendix the system of equations for m = —3, —1 will be presented, which were 
not investigated in Sect. IX. However, the method for their derivation is completely the 
same, but some new interesting consequences will appear. 

For the purpose, let us first rewrite the two sums in expression (111) for [p (z)] , 
putting in the first sum 2{n—l) — m and in the second term n — 4 = m. Then expression 
(111) acquires the form 

-, 2 4 i °° 

p'(z) =- + tt E G m+6 (m + 2) 2 (m + A) 2 z m - 

Z iD m=0,2,4,... 

-4 E G m+6 (m + 4)(m + 5)z m . (Al) 

m=-3,-2,-l,0,l,.. 

For m = —3 only the term from the second sum in (Al) will contribute. Putting also 
m = -3 in the R.H.S. of expression (118) for M(z)p 3 + N(z)p 2 + P{z)p + E(z) for the 
case of M = 4, N = 0, P(z) = -g 2 (z) and E(z) = -g 3 (z) (i.e. c^ = for all m , cj, 3) = 4 
and Cm = for m ^ 0),one has to take into account that terms with a " negative- valued" 
indice like c_{_ 2 „, c_[_ n) c_l_ 3n (n = 1,2....) are zero. 
The obtained equation for m = —3 is 

-8G 3 = 2(n + l)G n cfl n + c { l\ + c [ l\_ n G n + c [ % . (A2) 

For n — 1, when the first term on the R.H.S. is non-zero, the equation is 

-8G3 = I6G1 + c a i + c a lGi + c { % . (A3) 
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Since from expressions (142) and (136) from Section 7 c_g and c_{ can be found, from 
the above equation (A3) c_l can be expressed. 
For n = p > 1 equation (A2) is 

_8G 3 = <£} + £\_ v G p + c { % . (A4) 

Substracting the two equations (A3-A4), one can find the following expression for c_^_ p : 

■£Gp + 3(^) G p - (3^ + 21)Gp_ 2 



c_L. = ^r 16 + cli) = 16 — 



3-p /m l-^Ti^; ^ G 2/ ^ _ r« 



(A5) 



An expression for c_ 3 _ p can also be found from formulae (143) for k = p + 3 

16G 

_ GL 
r p+3 q 3 Lz p+ 



<&, = - r ~~ G : r ■ ( A6 ) 



From the equality of the two expressions, Gi can again be expressed as a convergent 
expression. Note also that the formulaes G 2p = G 3 (^ i )^^ and G 2p +i = G 3 (^ 1 ) p_1 from 
(151) satisfy the equality expression since then the denominators in (A5) and (A6) will 
be zero. This precludes the investigation of the system of equations for m = —3. 
For m — — 1, the general equation can be written as 

- - 48G 5 = 2(n + l)G n cfl n + (n + l)G n cfl n + Gjl\_ n + c? + cL°i , (A7) 

and for n = 1 and n = 2 the corresponding equations are 

-48G 5 = 8G1 + c2Gi + c[ 1} + cL°{ , (A8) 

-48G 5 = G 2 c { ll + ci 1} + c (0 l . (A9) 

The coefficient c_g can also be found from (143) for k = 3 



-,(i) 



160. 



U 3 — TtT^I 



(A10) 



G 2 

(1) , .(0) 



Substituting (A10) into (A9) gives an opportunity to express c\ + c_{ as 

c <« +c ™ = -48C 5 + 7 ^2g_. (A11) 



This expression, together with formulae (137) for c_ 2 and (143) for c_{, represented as 
c_[ = ^,can be substituted into the first equation (A8) to obtain the following quadratic 
equation for G\ 

(G 3 - ^±G 1 )(^F - 24G0 + 16G 2 G 4 = . (A12) 

Lr 2 Lr 2 
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2L = 


G2 Gk 1 

G3G4 Gfc-2 C3 




G 3~ g; Ui 



In a similar way, one can write down the equation for n = 3 

-48G 5 = 32G 3 + c^Gg + c? + c^} . (A13) 

Substituting c_\ and q + c_{ from (136) and (All), one can derive 

(A14) 

<_T3^4 ^fc-2 ^3 

where 

L = — G ^ . (A15) 

7^3(^3 — q|) + 2G2G4 

From the above two equations a relation, similar to (147) can be obtained 

G fe = /3G fc _ 2 = ^G 1 for fc = 2p 

= F>G X for fc = 2p + 1 , (A16) 

where 

0=^(21+1). (A17) 

Of course, in order to have an unique determination of Gfc,one has to require /3 — 7, where 
from (148) (3 = 2^ — ^-. This will result again in a quadratic equation for G\. 

Much more important and informative in the investigated case m = — 1 turns out to 
be the equation for a general n > 3 

-48G 5 = c { l\_ n G n + c? + cL°i . (A18) 

Let us remind that an expression for cl_[_ n can be written from (143) 

<£Li = ~r \r— ■ ( A19 ) 

Also, from (All) one has an expression for c\ + c_{. These two expressions can be 
substituted into equation (A18), which acquires the form 

G G 

G n (Gs — —Gi) = G2Gi(G n+ \ — — G n _i) . (A20) 

(j2 tj4 

Now it is interesting to note that using the formulaes 

G G 

G n +i = G2P+1 = Gz(—) p ; G n -i = G2p-i = G3{—) p " (A21) 

Cr2 Lt2 

from (151), it can easily be checked that the R.H. S. of (A20) is equal to zero for the case 
n = 2p. The other case n = 2p + 1 gives the same result. Therefore, from (A20) the 
following concise relation is obtained, expressing the proportionality of G\ and G2 with a 
coefficient of proportionality the ratio ^ 

G l = ^G 2 . (A22) 

(_T4 
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APPENDIX B: ADDITIONAL SYSTEM OF EQUA- 
TIONS FOR m=2k 



This Appendix will preclude the proof, started in Sect. IX that all the coefficient functions 
in the Loran function decomposition of the equation [p (2)] = 4p 3 — g2(z)p — g%{z) can 
be uniquely expressed, and especially those from the positive-power decomposition. 
For m = 2k > 0, the corresponding equation is 

(* + 2[(k + lf(k + 2) - 8(2k + 5)]G 2fc+6 = 2(n + l)G n cf (k+2) _ n + (n + 1) 2 G 2 J 2 % +1) _ 2n + 



+ (n + l)G n c 2{k+1 y n + 2(n+l) G 2n c 2 ^ k+l) _ 2n + {n + l) G 3n c 2k _ 3n + c 2{k+1) +c 2k _ n G n + c 2k . 

(Bl) 
Additionally fixing the value of n — k + 1, one can obtain from (Bl) 

(A;+2)[(A;+l) 2 (A;+2)-8(2A;+5)]GW 6 = 12(k+2) 2 G 2{k+1) +c^ k+1) +G k+1 c k 1 ! 1 +c 2 ° k ) . (B2) 
For n = 2(k + 2) the equation is 

(A;+2)[(A;+l) 2 (A;+2)-8(2A;+5)]GW 6 = %{2k+^)G 2{k+2) +c { 2 % +1) +c^\G 2{k+2) +cf k ] . (B3) 

Substracting the two equations, one can express c k _ x as 

4-i = ?r-[-12(A; + 2fG 2{k+1) + 8(2A; + 5)G 2(fe+2) + c W A G 2{k+2) ] . (B4) 

Since fc > 1, from this formulae it is clear that all the Loran coefficients c m in the positive- 
power decomposition can be expressed, including the coefficient c , through which the 
coefficients from the negative-power decomposition in Sect. IX were expressed. Also, 
from (B4) cLjj, +1 -, can be expressed (by performing the indice change k — 1 — > 2 [k + 1). If 

c 2(i+i) * s substituted back into equation (B3), one can express also the even positive-power 

coefficients c 2k as 

eg? = (k + 2)[(k + l) 2 (k + 2) - 16(fc + 3)]G 2fc+6 - [cL 1 ! + 8(2A; + 5)]G 2(fc+2) + 

--4 

Now let us write down equation (Bl) for another possible value of n — 2(k + 1 

2(fc+l) 



f^ [12(2fc + 5) 2 G 4 ( fe+2 ) - (32A; + 88 + c^)G 2{2k+5) ] . (B5) 

<J"2(fc+2) 



(A;+2)[(A;+l) 2 (A;+2)-8(2A;+5)]G' 2fe+6 = 4(2k+3)G 2{k+1) +c^ k+1) +cf k ) +c W 2 G 2{k+1) . (B6) 
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Combining (B3) and (B6), c_ 2 can be expressed as 

c (}) 2 = — — [-4(2* + 3)G 2(fc+1) + 8(2A; + 5)G 2(fe+2) + c^ A G 2{k+2) ] . (B7) 

t*2(fc+l) 

The coefficient c_\ can easily be calculated from (136) and (142) to be 



16 
& = - 



3-rr\G 2 Gk — G^Gk-2) — 2QG A G k - 2 



G2Gk — G^Gk-2 



(B8) 



On the other hand, it is important to observe that c 2 can be calculated independently 
also from equations (100) and (142) 

(B9) 



(i) _ ^\G 2 Gk — G^Gk-2) — I6G4GA 



^T^k ~ ( - T 4 < - T fc-2 

Note also that from relations (151) for G*2p and (?2p+i it follows 

G G 

-J^-l = _i for k = 2p and k = 2p + 1 (BIO) 

Gk G4 

or written in another way - G 2 G k — G^Gk-2 = 0. 

Setting up equal the two expressions (B7) and (B9) for c_ 2 , cancelling the equal 
denominators and subsequently taking into account (BIO), one can obtain the following 
concise recurrent relation 

G k = a(k)G k . 2 , (Bll) 

where a(k) denotes 

a(A;)=20^^. (B12) 

G2(k+1) 

Continuing further the recurrent relation (Bll), one can derive 

G k = a(k)a(k - 2) a(k - (k - 3))G k -( k -i) = 

tt G 2 ( k+ 2) G 2k G 2 (k-2) G2.5 /t3iq\ 

= coeff.- ....__Gi . (B13) 

( - r 2(fe+l) ( - r 2(fc-l) ( - r 2(fc-3) ^2.4 

If k = 2p, the numerical coefficient in (B13) will be 20 p . 

A similar relation can be obtained by fixing n = 2(k + 1). Then the corresponding 
equation is 

(k + 2)[(k + l) 2 (k + 2) - 8(2k + 5)]G 2k +6 = 

= A(2k + 3)G 2(fc+1) + <£{G 2(Jfc+ i) + c% +1) + cg> . (B14) 
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Substracting this equation from (B16) for n = 2(k + 1) and taking into account expression 
(B4) for c k _ v one can obtain 



Gk-l 



cL 1 i + 8(2A; + 5)]G 2(fc+2) 



= G 2{k+1) [<£> + A{2k + 3) + 12(fc + 2) 2 (^ + i g fe-i) j (B15) 

where c_4 is given by (B8) and c_{ by (B9). Substituting the above expressions into 

(B15) and again taking into account that G 2 Gfc — G^Gk-2 = 0, one derives the following 

recurrent relation 

Gk-i 

Gk+i 
If this relation is substituted into (Bll), then it can be derived that 



G 2{ k+i) = 20^ . (B16) 



G k Gk+i — 20Gk-iGk-2Gk+2 ■ (B17) 

This equality is valid for k > 3. For k — 3, 5, 7, 9 the above relation may be written as 

GrgCrg = 20Gltr5 ^5^12 = Gr^Grj , (B18) 



^7^16 — ^5^9 ^9^20 — G7G11 . (B19) 

Since on the L.H.S. of (B18) and (B19) Gg, Gu, G±e and G20 are zero, the R. H. S. should 
also be zero. If G\ 7^ 0, G3 7^ 0, the R. H. S. of the first pair of equations (B18) equals to 
zero if G5 = G-j = 0. But since G 5 and Gj appear also in the R.H.S. of the second pair 
of equations (B19), the R.H.S. will be zero and therefore Gg and Gn may be different 
from zero. The treatment of the subsequent equations is analogous. That is why one may 
conclude that a pair of even sums G21+1, G21+3 (I > 2) is zero, but the next pair G21+5, 
G21+7 may be different from zero. 

The last fixing of the value of n = ^ for the case m = 2k gives the equation 



(A; + 2)[(A; + l) 2 (A; + 2)-8(2A; + 5)]G' 2fe+6 = 4(l + ^) 3 G' 2fc + c« +4 ) G^+c 2 ° fc ) . (B20) 

O 3 3 



Substracting from (B20) equation (B6) for n = 2{k + 1) and setting up y = p, one can 



,(1) 



\3/n. Jl)/ 



12(p + 1) + c^j G 3p+2 - 4(1 + pYG 3p - c^G p = . (B21) 

Similarly, substracting from (B20) equation (B3) for n = 2(k + 2), one obtains 

4(1 + p) 3 G 3p + <MG P - 8(3p + 5)G 3p+ 4 - c { l\G^ +A = . (B22) 
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From the two equations it follows 



.(i) 



(i) 



12(p + 1) + c^J G 3p+2 = [8(3p + 5) + cLlJ G 3p+4 . (B23) 

Again taking into account (B8) for c_ 2 and (B8) for c_ A for value of k — 3p, one can 
obtain 

G 3p+4 = 20^^G 3p . (B24) 

^3 P -2 

However, in view of the relations (B18-19) and the consequences from them, the last 
relation will make sense only when each of the indices 3p + 4, 3p + 2, 3p — 2, 3p equals one 
of the indices 2/ + 5, 11 + 7 (I > 2) and then the relation (B24) will be nonzero. 

APPENDIX C: ADDITIONAL SYSTEM OF EQUA- 
TIONS FOR m=2k+l AND m=-k 



For m = Ik + 1 the corresponding equation is 

-8(2* + 5)(fc + 3)G 2fc+7 = 4J +5 + 2(n + l^cg^ + (« + l) 2 ^^.^ 
+(n + l)G n cS + 3-n + 2(n + l) 2 G 2 J$ +3 _ 2n + (n + l) 3 G 3 „4 3 fc ) +1 _ 3 n+ 

, Ji) , Ji) r + r (0) cm 

"r c 2fc+3 "r ^k+l-n^n ~T~ c 2£+l ■ l^-U 

The important conclusion, which can be made from this equation is the following: if c 2k+3 

is calculated from (B4) for value of k — 1 = 2k + 3, then the odd number coefficients c 2k+1 
can also be found! Remember also that in Appendix B only the even number coefficients 
c 2k were found (form. B5). In order to express c 2fe+1 , it is enough to set up n = 2k + 3 

in (CI), when from all the coefficients c m only the second term on the R.H.S. will be 
non-zero. Then 

41+1 = -4+3 - c { l\G 2k+5 - 8(fc + 3)(2k + 5)GWr - 16(fc + 3)G 2k+5 . (C2) 

For another value of n — ^+j^, equation (CI) acquires the following form 

-8((fc + 3)(2k + h)G 2k+7 = 3(2k + 5) 2 G 2fe+3 + cg+a + 4^^ + cf k ] +1 . (C3) 



2 



But since c 2fc _i and Ga+3 have to be integer numbers, this will be possible if for example 
2k — 1 = 2p. For this value of k, one can express c 2 +2 from (C3) 



C 2p+2 



"" - -4p+4 - 4 1)g p+2 - l2 (P + 3) 2 G 2p+4 - 8(p + 3)(2p + 7)G 2p +s . (C4) 
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However, c 2p+4 can be expressed also from equation (B2) for values of m — 2k, n = k + 1 

and k = p + 2 



„(o) _ _ (1) 

^2p+2 ~~ C 2p+4 



J p 1) G p+2 -12(p + 3) 2 G 2p+i + 



+(P + 3) [(p + 3)(p + 2) 2 - 8(2p + 7)] G 2p+8 . (C5) 

From the two equations (C4-C5), one easily obtains 

lQ(p + 3)(2p + 7)G 2p+8 = 0. (C6) 

Since the coefficient in front of G 2 p+s is a positive one, (C6) will be fulfilled if 

G 2p +s — G 2 k+e = . (C7) 

The last means that the even-number sums G 6 , Gg, Gio, are zero! 

Again fixing the value of n = 2k + 3, one derives from (CI) the equation 

-8(fc + 3)(2fc + 5)G 2fe+7 = [8(k + 2) + c^ 2 ]G 2k+3 + c^ +3 + c£ } +1 . (C8) 

Combining this equation with (C3), setting up 2k — 1 = 2p, one can express c p 

cf = _W)-^)Q^ , (09) 

where 

A(p) = 3(2p+l) 2 + 26(2p+l) + 59 . (CIO) 

Comparing this expression with formulae (B4) for c p and taking into account (B8) and 
(B9), one derives the following relation 

n 1 ^2(p+2) ,„ s 

G 2(p+3) = 7^-^ • l^ 11 ) 

ZU Or 2 (p-2) 

Expressing by means of (B15) G 2 ^ p+2 ) and G 2 ( p -i) and substituting into (Cll), one can 
obtain the relation also in another form 



1 G p + 2 G p _a 

20 Cttj Gr r , 



^2(p+3) — 7^{— — ~Fi > X^ 11 ) 



and from the two expressions one can obtain also the ratio 2(p+2) without any numerical 

L Cr 2 ( p _2) J 

coefficients. Note also that (C11-C12) refer to non-zero even numbers of G m since we 
have 2p = 2k — 1, and the relations should be written in respect to k and not p. 
The corresponding equation is 

-8(k + 3)(2fc + 5)G 2k+7 = ^-(2k + 4) 3 G 2fc+1 + c ( 2k \ 3 + G^l]^ + c^ +l . (C13) 

Z I 3 3 
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Setting up 2k + 1 = 3p and keeping in mind that c 2p and cL +2 can be found from (B4), 

(o) 
one can express c 3p as 

4p = ~^(3p + 4)(3p + 5)G 3p+6 - 4(p + l) 3 G 3p - cg +2 - G p c% . (C14) 

The equation for the last case of m — —k (k > 3, k ^ 6) is 

cf\ + 2(n + l)C n cL 3 l +4 _ n + (n + l) 2 G 2 „c£ 3 i +2 _ 2n + (n + l)G„c<_ 3 J + 2-n+ 



+2(n + l) 2 G 2n c^ +2 _ 2n + (n + l) 3 G 3 „c {3 L 3n + c« +2 + c ( _l n G n + c^ = . (C15) 

Since it has been shown already how all the Loran coefficient functions can be expressed 
and the treatment of this equation is completely analogous to the preceeding ones, equa- 
tion (C15) shall not be considered. 



1 APPENDIX D: COEFFICIENT FUNCTIONS 

N U N 2 ,N 3 and N 4 DEPENDING 
ON THE "BAR" VARIABLES 



The coefficient functions N\, N2, iV 3 and N4 in the cubic algebraic equation (106) for T in 
Section VII are the following 

N, = 2p\Z 2 - 2Z%(^ + ^p 2 + p 3 ) - Z%(0 2 p, + Op 2 + 2p 3 ) + 

+Z«(0 2 p 1 + Op 2 + 2p 3 )(^l + °EL +p 3 ) (Dl) 

N 2 = 8Zpj - Z%(p 2 + Op,) + [Z\p 2 + Op,) - lGZ 3 ^] (^ + ^p 2 + p 3 ) + 



+(0 2 Pi + Op 2 + 



15 , 23^ 5/ 2 _ 0_ _, 
— Z 3 Pl + —Z\— Vx + -p 2 +p 3 ) 



(D2) 



N 3 = (^U + % 2 + p 3 ) [^Z\0 2 p x + Op 2 + 2p 3 ) + 8Z 3 (p 2 + Oft) - 2Ap,Z 2 ] 



I ' ' 2" 

+8p 2 - 4Zp 1 (p 2 + OpO - 8Z 2 p 1 (0 2 p 1 + Op 2 + 2p 3 ) (D3) 

iV 4 = -4ft fo + Op x ) + 6^(0^ + Op 2 + 2p 3 ) + 
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+ [-18^(0^ + 0p 2 + 2p 3 ) + 12Z 2 (p 2 + Op,)] (^ + ^p 2 + p 3 ) . (D4) 
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